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GROWTH AND DIFFERENTIATION IN LIVING PLANT 
MERISTEMS 


By Epmunp W. SINNOTT 
DEPARTMENT OF Botany, COLUMBIA UNIVERSITY 


Communicated January 12, 1939 


Most of our knowledge of the cellular changes taking place during growth 
has been derived from a study of killed and sectioned material. Cell 
division in living plant tissues has been observed in a considerable number 
of cases, but chiefly in hairs, protonemata and similar structures which 
are only one or two cells thick. It is desirable to learn what takes place, 
in terms of changes in cell number, size and shape, in a typical living plant 
meristem as growth proceeds. Events cannot readily be observed here 
because most meristems are so massive as to be opaque. Furthermore, 
the terminal growing point of a stem is buried in surrounding structures, 
and that of a root is usually enclosed in a relatively opaque calyptra. 

These difficulties may be overcome by the use of appropriate material 
and technique. Seedling roots growing from very small seeds are so 
delicate, often only six or eight cells thick, that light passes through them 
readily. In most cases, however, direct observation of the actual meri- 
stem is prevented by the root cap; but in some of the small-seeded grasses 
the cap is reduced to a group of only a few rather loose cells and the meri- 
stem is thus nearly naked. Species of Agrostis, Phleum, Poa, Sporobolus 
and related genera provide particularly favorable material. 

It has been found possible to grow these very delicate roots under 
essentially normal conditions by germinating the seeds on strips of moist 
lens paper on microscopic slides, placed in covered staining dishes. About 
a centimeter of tap water is left in the bottom of the dish, which keeps 
the air humid and the paper moist. After about four days at room tem- 
perature the young roots are usually more than a centimeter long and are 
ready for study. The slide is now placed on the stage of a microscope 
and the root covered with a few drops of water. The meristem may be 
examined directly under a water immersion lens, but more satisfactory 
results have been obtained by placing a cover glass over the root, sup- 
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ported by glass chips on either side. The root may now be studied under 
any immersion objective, but the author has found water more satisfactory 
than oil for contact, since the latter tends to spread over the edges of the 
cover. Care must be used to avoid any pressure on the root when the 
cover glass is applied or removed, and to keep the root continually moist, 
for otherwise growth will stop immediately. By means of a camera 
lucida, drawings may be made of the cell outlines of the surface layer of 
cells, the second layer and sometimes the third layer. The slide should 
be returned to the damp chamber as soon as possible. 

Drawings may be made at intervals of from one to several hours. If 
a marker, such as a cell of unusual shape or size, is located in the first 
drawing, then all the cells may be numbered with reference to it, and 
located again in the second and later drawings. In this way, all the new 
cell walls which have appeared through division in the interval may be 
recognized. Changes in cell size, both before and after division ceases, 
may also be measured. The actual cell lineage for a long series of cells 
from the apex of the meristem to the point where elongation has ceased 
may thus be determined. The nuclei are to be seen only vaguely and 
details of mitosis are not visible. 

These roots are very sensitive to gravity and will not continue to grow 
normally more than a few minutes on the horizontal stage. It is more 
satisfactory in some respects to keep them vertical, with the microscope 
in a horizontal position, but this makes difficulty in the maintenance of 
liquid contact between lens and cover, and in the use of the camera lucida. 
Photographic records may be made, but unless conditions are very favor- 
able, these have been found less satisfactory than drawings. The chief 
difficulty thus far encountered is the tendency for the growing tip to twist 
spirally and thus sometimes to prevent continuous observation of a given 
row of cells. 

The approximate rates of cell division and of cell elongation may readily 
be measured. The rate of elongation of the whole root may also be ac- 
curately determined by marking the location of a given point under the 
camera lucida at constant time intervals. The characteristic rhythmic 
nutation of the root tip may be observed and its period and amplitude 
measured. 

The method here described makes possible a fresh approach to the study 
of a number of other problems of growth and differentiation in a typical 
terminal meristem. Preliminary results in three of these problems, as 
studied in the genera Poa, Phleum, Agrostis, Chloris and Sporobolus, are 
presented briefly here. 

1. Factors Determining the Shape and Position of New Cell Walls.— 
The cells in all layers are in regular longitudinal rows, so that the new 
walls, with rare exceptions, are transverse or essentially so. In equa- 
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tional divisions the walls are straight, but in non-equational ones, as in the 
divisions giving rise to root-hair initials in some genera, they are com- 
monly curved from the first, with the convex side toward the larger cell. 
A new wall is never laid down directly opposite a previously formed cross 
wall of any neighboring cell. If the natural position for the new cross 
wall would tend to bring it into such a relation, this will be avoided by a 
bending of the new wall so that it joins the side wall at a point at one side 
or the other of the place where the cross wall of the neighboring cell is 
attached. A new cross wall tends to join the side wall halfway between 
the points of attachment of the two nearest cross walls, either of this cell 
or a neighboring cell. This is the position which would be assumed if the 
walls were liquid films. In many cases, however, this position would result 
in daughter cells very unequal in size, and in such cases the new wall is 
pushed considerably nearer the point of attachment of one of the cross 
walls. The point where the new cross wall joins the side wall thus seems 
to be determined by an equilibrium between the tendency to divide the 
mother cell equally (or into two parts of definite proportions) and the 
tendency to bisect the free side-wall space. 

2. ‘Sliding Growth.’’—Since the point where each new cell wall is laid 
down can thus be definitely located with reference to the walls of all ad- 
jacent cells, and since this relationship can be followed throughout the 
history of all these cells, it is possible to determine with some certainty 
whether the cells slide along each other during the process of growth. This 
may best be studied in the relationships between the cells of the surface 
layer. All the evidence indicates that sliding growth does not occur. The 
relative positions of the transverse walls of adjacent cells are maintained 
essentially unchanged from the time they are laid down until the cells are 
many times their original size and growth stops. There is a definite 
gradient in cell expansion, beginning just before division ceases, and 
gradually increasing until maximum growth is attained. In this all the 
cells participate almost equally, though the root hair initials and certain 
stunted cells grow somewhat less rapidly than the rest. The relative 
position of their walls thus changes a little, but this is clearly due to a 
greater growth in some parts of the wall than in others and not to sliding 
growth. No case has been observed where the end wall of one cell has 
passed by the end wall of an adjacent cell during growth. The relative 
wall positions between the surface layer of cells and the layer immediately 
below have also been determined and here, too, the positions established 
when the cell walls are laid down persist to maturity. 

3. The Differentiation of Root Hairs—In Phleum, Poa and Agrostis, 
the cells which are to produce root hairs (the “‘trichoblasts” of Leavitt) 
are definitely set apart at the last cell division. This is usually either the 
third or the fourth division from the tip of the root. At this division the 
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apical (distal) member of the pair of daughter cells is always smaller than 
the basal. It has denser protoplasmic contents and its walls tend to bulge 
outward slightly. At some distance back from the tip, when elongation 
has ceased, this cell produces a root hair. The other cell, which rarely 
may divide again, never produces a hair. In these genera, therefore, 
differentiation for this character occurs very early and the potencies of 
the cells are sharply limited almost from the beginning. In Chloris and Spo- 
robolus, on the other hand, trichoblasts ordinarily do not appear, all the cells 
being essentially equal in size from the beginning and all being capable of 
producing root hairs. Some of the cells stain differently from others, and 
a differentiation of root hair cells may occasionally be observed before the 
hairs are formed. These plants evidently provide good material for a 
study of the factors controlling the differentiation of root hairs. 

Summary.—Materials and methods are described, by the use of which 
it is possible to observe and measure the multiplication and growth of 
cells in living root meristems. This technique has been applied to the 
problems of the location of new cell walls, of ‘‘sliding growth,” and of the 
differentiation of cells which are to form root hairs. 

1R.G. Leavitt, Proc. Boston Soc. Nat. Hist., 31, 273-313 (1904). 


QUANTITATIVE ANALYSIS OF THE INTERACTION OF 
INDIVIDUALS 


By Extot D. CHAPPLE 


DEPARTMENT OF ANTHROPOLOGY AND DEPARTMENT OF INDUSTRIAL RESEARCH, HARVARD 
UNIVERSITY 


Communicated January 13, 1939 


Up to the present time, no quantitative studies of human interaction 
have been undertaken, with the exception of observations by Dorothy 
Thomas and her co-workers on the behavior of children.1 These studies 
were designed to develop an “‘index of personality” based on the percentage 
of total iime spent by a child handling objects, interacting with people or 
playing alone. Only percentage figures were secured for interaction, and 
the authors were more interested in developing criteria for determining 
the reliability of observers. 

This paper represents a preliminary account of quantitative results 
secured through the use of a crude recording apparatus, now being super- 
seded by an accurate instrument. Although the investigation was pri- 
marily exploratory, it is believed that a discussion of the results obtained 
by use of the old instrument will be of interest as an indication of the 
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possibilities of research in this field. It must be emphasized that the 
conclusions herein presented are entirely tentative. They must be tested 
by the use of accurate apparatus under laboratory conditions (using a 
one-way screen, etc.). This work is now in progress. 

The present studies resulted from the development of operational 
methods for the description of the relations of individuals. These opera- 
tions were evolved in the course of the analysis of material collected in 
field studies of communities. From the consideration of this material, a 
working definition of the class of phenomena was derived, which may be 
stated briefly as follows: Individuals are considered to be in interaction 
if the action of one individual is followed by the action of another individual. 
The full definition and a description of the analytic procedures are given 
elsewhere.? Here it is sufficient to point out that such an action is a 
manifest (hence observable) phenomenon, and may be made up of words, 
gestures, in general, of overt muscular activities. No distinction is made 
between kinds of actions. 

In order to obtain a record of a sequence of actions manifested by 
individuals, it is necessary to use some kind of a time-recording apparatus. 
Accordingly, a simple device was improvised. A large wheel was fitted 
to the rubber roll of an old noiseless typewriter. A small electric motor 
drove this wheel by friction at a uniform rate of speed (15 inches to the 
minute). A roll of adding-machine paper, mounted on a brass frame, 
was fed through the roll to a take-up driven by another small motor. The 
take-up motor rewound the paper on a wooden roll and was arranged so 
that no pull was exerted on the paper coming through the typewriter. 
Thus, any influence on the speed due to the changing size of the upper 
roll was eliminated. An observer, seated at the typewriter, struck a 
designated key when the first individual acted, and when this action was 
ended by the action of the second individual, another key was struck. 
By the alternation of the two letters on the moving tape, a continuous 
record of a conversation could be secured. The length of time of each 
action was obtained by measuring from the top of the first letter to the 
top of the second. We thus obtained a series of durations of the actions 
of two individuals. 

Observations were made only on two individuals at a time, due to the 
limitations of the machine. Other difficulties in recording, now elimi- 
nated in the new apparatus mentioned, were due to the inability to record 
those instances when two individuals talked at the same time. (In this 
case, individual B was recorded as ending A’s actions, and then almost 
immediately A comes in again.) It was also impossible to mark where 
one individual fell silent in the midst of an action and began again when 
the other individual was not observed to manifest an action. Because 
of the lack of proper facilities for observation the observer had to be seated 
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in the same room, with the subjects at a distance of eight to ten feet away. 
The effect of this situation was to a large degree minimized, due to the 
fact that the subjects were colleagues who had known each other for several 
years, and the conversations were those which would ordinarily have 
taken place. The same observer was used in all observations described 
in this paper. For this reason, the error due to the observer may be re- 
garded as constant, and is neglected in the rough approximations here 
described. We shall consider two series of eight observations, and one of 
five. One individual was a member of all three pairs. Due to difficulties 
in getting the subjects together, only a few observations were obtained 
on pairs made up of the other three individuals, and we shall therefore dis- 
regard those figures. 

In each observation, lasting from thirty to fifty minutes, the two indi- 
viduals manifested a total number of actions ranging from 320 actions (155 
interactions) to 840 actions (420 interactions). The duratidns of any 
single action varied from 0.2 seconds (the minimum that could be recorded 
on the machine) to extreme values of well over a minute. Frequency dis- 
tributions showed a characteristic asymmetry, the curve being definitely 
J-shaped. When one-second intervals were used, the mode fell in the 
interval between 0.2 and 1.1 seconds, containing about one-third of all 
values. The means of different individuals differ from observation to 
observation, the lowest 2.83, the highest 7.81, but approximately two- 
thirds of the values were found between the mean and 0.2, the other one- 
third in the higher values. In the present paper, we shall concern our- 
selves not with the probability that any single unit of action will have a 
given duration, but rather with the way in which the sequence of dura- 
tions is arranged. 

In order to describe the order of occurrence of the action-durations of 
two individuals, we obtained a kind’ of internal average by summing each 
five units of action. For reasons that need not be discussed here, attempts 
to use running averages were discarded. This summing of five should 
average out errors due to the observer; it is sufficiently small to be sensitive 
to fluctuations, yet at the same time affords the advantages of an average 
while enabling us to obtain the sum of the durations of any given series 
of actions. By comparison of the variation of the values of both indi- 
viduals, A and B, we could begin to estimate the effect of one on the other. 
A visual representation was afforded by plotting the sequence of sums of 
fives for each individual on a graph, the ordinate representing the cumula- 
tive number of actions and the abscissa the cumulative time (sum of five) 
of the actions of each person. The same procedure was used in treating 
the units of interaction, summing by fives and plotting the results cumu- 
latively. 

Upon inspection of the curves for the separate individuals, it was seen 
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that no single empirical function described all the curves in a single series 
of observations. We therefore provisionally regarded the curves as made 
up of a series of straight lines, each describing that period of action during 
which an individual’s rate was approximately constant. By selecting 
points on the curve in which the variation was a minimum, we calculated 
the slopes between the points and thus obtained a sequence of slopes for 
each individual. The rate of any individual for a given slope is defined 
by 2t/n. 

It should be pointed out that in the present study we have to regard 
these slopes as making up a sequence of rates of silence for each individual. 
This way of looking at the facts derives from the operations used to obtain 
the discrete durations. With the particular apparatus used, a typewriter 
with a moving tape, we are unable to obtain a measure of the duration of 
an action of one individual uninterrupted by the action of the other 
individual. Let us suppose that one individual is observed to act con- 
tinuously without pause for twenty seconds. If the other person acts 
by nodding his head or uttering a monosyllable, this single action would 
be recorded as made up of three or four shorter actions. On the other 
hand, if one individual remains silent while the other individual continues 
to talk, the silent individual obviously cannot be interrupted. His 
silence will continue until he manifests an action. In other words, with 
this machine we could not distinguish the pauses an individual makes in 
an action during which the other person acts from those instances when 
both were talking at the same time. If the reader wishes to regard these 
figures as representing rates of action, he merely needs to shift the figures 
about. It should be pointed out, however, that this way of regarding 
the figures is not justified by the operations, and the results obtained 
(now to be described) reinforce this view. 

In table 1, we present the slopes (mean rate), the sum of the durations, 
as well as the durations of each slope in chronological time for two con- 
versations from different series. The table is arranged in order of the 
occurrence of these slopes in the sequence of units of interaction. From 
this table, it may be seen how one individual often maintains a slope 
while the other changes several times. There is no simple relation be- 
tween the durations and the values of the slopes; nor is there any evident 
correlation between the slopes of one individual and the other. 

In order to obtain an approximate description of the interaction of 
individuals represented by these slopes, we shall have to consider three 
questions: (1) what determines the value of the slope of an individual; 
(2) what determines the duration of any single slope of an individual; 
and (3) what is the nature of the dependence of the individuals upon each 
other; that is, is this relationship expressed in the associated slopes and 
their durations? 


) 

if 
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Let us first take up the question of the relationship of these slopes to 
one another. It is evident that we cannot express a single slope of indi- 
vidual A as a function of the slope of individual B, since we do not find 

TABLE 1 
SLOPES AND DURATIONS OF SLOPES IN Two EVENTS 


Observation No. 37 


NUMBER OF suM oF EDC’s DURATIONS Recnng GB's DURATION 
ACTIONS IN EDC ACTION PER OF EDC’s GB's ACTION PER or GB's 
SLOPE SLOPE SLOPE SLOPE SLOPE SLOPE SLOPE 

0-10 6.11 61.1 93.5 3.65 54.7 131.4 
10-15 2.49 62.3 437.6 
15-25 29.01 290.1 318.4 
25-35 6.17 339.1 761.3 
35-55 9.69 193.9 303.8 
55-70 3.68 55.2 157.7 at 
70-75 21.66 108.4 132.1 
75-80 7.73 773.0 1345.3 = = a 
80-115 re 2.78 97.2 393.0 

115-160 4.75 213.7 566.8 
160-175 14.76 221.4 291.1 


Observation No. 29 


SUM OF TIME 
NUMBER OF sum oF EDC’s DURATIONS CMA’s DURATION 


ACTIONS IN EDC ACTION PER or EDC's CMA’s_ ACTION PER or CMA’s 
SLOPE SLOPE SLOPE SLOPE SLOPE SLOPE SLOPE 
0-30 3.15 94.5 166.6 2.57 128.4 257.9 
30-50 1.96 78.3 229.6 
50-70 3.47 243.6 414.7 
70-85 3.71 55.7 90.2 de 
85-105 1.06 21.1 83.8 
105-120 3.65 292.2 559.1 
120-150 2.22 66.7 188.9 
150-155 10.84 54.2 72.5 
155-180 2.50 62.5 143.7 
180-185 9.06 90.6 113.9 
185-190 2.84 298.7 612.1 
190-240 2.06 103.1 282.4 
240-290 3.04 151.8 267.4 


Note: The broken line, __, indicates that the above slope continues. It should be 
noted that EDC is characterized by a pattern made up of a slow slope followed by a fast 
one, a simple alternation, while GB has wide variations due to two very slow slopes, 
and CMA exhibits a fast-slow pattern but with intermediate pairs of values having a 
different ratio than the first and last pairs exhibited. 


one slope of individual A associated solely with a single slope of individual 
B. In fact, such a situation is very rare in this material (five times in 
333 slopes manifested); while one individual maintains a constant rate, 
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the other individual may manifest as many as five different slopes. If 
we hope to describe the duration of a single slope of individual A as some 
function of the slopes that B manifests during the time that A is constant 
we shall have to express the relationship as one between the durations. 
By estimating the ratio of the durations of A’s slope to the total time 
(the sum of the durations of the slopes of both A and B), we can then see 
whether this varies as a function of the value of the slope. 

When log slope is plotted as a function of this ratio for all the slopes 
manifested by one individual in interaction with another individual, it 
is found that the points representing the relationship make up a band. 
The band has a constant width and may be enclosed by two parallel lines 


50 
100 (Tj/ Tat Tp) 
FIGURE 1 


fitted to the points making up the border. Figure 1 represents the band of 
GB with EDC. 

When the bands for the three pairs discussed here were plotted, it was 
found in the equation 


log Sa = k(ta/ta + ts) + constant 


where S equals the slope, and A and B are the two individuals, and ¢ the 
sum of the durations, that the constant describing the slope of the band 
was the same for individual EDC in the three pairs in which he interacted, 
differing in slope from the bands of GB, CMA and FLWR, who also differed 
from each other. Unfortunately we were unable to obtain a sufficient 
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series for the other individuals with each other, but it is an obvious possi- 
bility that the constant & in this relationship is invariant for each individual. 
This similarity in the value indicates why these slopes should be regarded 
as representing silences. Since the same slope is found for the silences 
and not for the individual’s actions in three separate series, it may be 
assumed either that the individual adapts the durations of his actions to 
the other person and does not adapt his silences, or that the recording 
device’s limitations are responsible. No solution is possible, of course, 
until a more accurate apparatus is employed. 

What precisely does the band function tell us? In the first place, it 
is evident that the function imposes limits upon the duration of any single 
slope. If individual A manifests a slope of a particular value, and B 
manifests a slope which does not bring the value of the ratio within the 
limits of the band at the given value of log slope for A, then individual A 
will maintain this slope, until other slopes manifested by B bring the 
ratio within the band limits. Thus, the duration of a slope at a given 
value is determined (within certain limits) by the proportion it makes 
with the total time. For the different individuals with whom EDC inter- 
acts, however, the band for the function has a different width, and it is 
suggested that the width of the band and the distribution of points within 
it is a measure of the degree of adaptation of the individuals to one another. 

Not only do individuals seem to possess constant slopes with bands 
having definable limits; in one case, an individual, FLWR, exhibits a 
marked discontinuity. At a value of log 0.4800, the whole band shifts 
towards the X-axis and changes its slope, although the width of the band 
remains constant. This indicates, for this individual at least, that when 
he slows down to rates of silence above 3.02 seconds, the duration of his 
slope forms a greater proportion of the total time, and he has, therefore, 
differing rates of adjustment in interaction. 

It should be pointed out that the limits of these bands cannot properly 
be regarded as defining the probability that repeated measurements of 
log slope as a function of the ratio will fall within the band instead of 
outside it. We shall come upon such a function a little later in this paper. 
In this case, however, the situation is more complex. If a band encloses 
the values of the function for each conversation, it will be observed that 
the width varies from day to day, although the slope remains constant. 
In some observations, all the points fall on a straight line; in others, points 
forming both borders of the total band are included in a single observation. 
The question of probability, then, involves the distribution of the varia- 
tion for single observations. The limits of the band for all observations 
on a given pair define this probability. The view that the probability 
distribution is made up of bands for single observations is supported by 
the fact that the limits of the band do not act as limits which rigidly 
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control the change of slope. That is, when individual A changes to a 
certain slope, the first slope of individual B may bring the ratio within 
the band limits. That individual A does not shift depends upon two 
factors, (1) the existence of minimum and possibly maximum durations 
for slopes, the former giving a kind of threshold effect, the conditions for 
which have not yet been investigated; and (2) the limits of variability 
defined by the width of the band in a single observation. In this latter 
case, the variation in width of the bands is an expression of the variation 
in adjustment of the slopes of one individual to the slopes of the other, 
and hence depends upon the variability of the constituent individuals. 

If we are to investigate this variability further, we must develop some 
way of dealing with the individual values of the slopes, regarding the 
preliminary definition of the problem of their durations as sufficient for 
our present purposes. In one-half of all observations, 13 out of 23 for 
EDC of which seven were with CMA, and a smaller number for the other 
individuals, the sequence of slopes forms a pattern of alternating fast and 
slow slopes. In the majority of these cases, there is a constant relation- 
ship between the fast and slow slopes. That is, if each sequence is broken 
up into a series of pairs, and each slow slope in a pair divided by its ac- 
companying fast slope, the ratio or its first difference is constant for all 
pairs in the observation. The constant varies from observation to ob- 
servation. If we assume provisionally that this alternation of fast and 
slow slopes is characteristic of individuals in interaction, we must then 
try to show why deviations from this pattern take place. We shall not 
attempt to determine the conditions defining the ratios of the slopes nor 
the variation in their absolute values when the ratios are constant. 

If a consistent explanation is to be introduced to explain the deviations, 
we must find out whether the variability in the values of slopes differs 
for each individual or whether all the individuals under consideration 
vary in the same way. In order to obtain a measure of this variability, 
the dispersion of the slopes was expressed as a function of their mean value. 
If log o slope is plotted as a function of log Mean slope, where ¢ = 
—V/=x?/n — 1, the measurements for all four individuals fall within the 
limits of a narrow band defining the probability that measurements of o 
as a function of the Mean will fall within the band instead of outside it.* 
The slope is greater than one and the origin is on the abscissa. Hence 
log ¢ = k log M + Constant. It should be pointed out that a plot of o 
of the raw figures against the Mean of the raw figures gives us a similar 
band with the slope approaching 1. It is, of course, premature to suppose 
that this function will be characteristic of all individuals on the basis of 
such a small sample, but it is of interest to point out that similar results 
have been obtained by Crozier and his co-workers for a wide variety of 
biological material.‘ 
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Since we know that these four individuals vary in exactly the same way, 
we do not have to regard the occurrence of extreme slopes for any single 
individual as due to a special property of that individual. Rather this 
variation or the causes thereof would seem to be common to all the indi- 
viduals studied. If we assume that every individual alternates between 
fast and slow rates, then deviations from this pattern must have a con- 
sistent explanation. These deviations are uniformly very slow rates 
occurring either in a position where a slow rate was expected, but not 
one so large (on the theory that in any conversation, if the ratio of each fast 
slope to its accompanying slow slope is calculated, the constant derived is 
constant for all such pairs in a single conversation) or else such deviations 
occur in place of a fast slope, thus disturbing the regularity of the pattern. 
If the rates and durations of the other individual are then examined to 
see whether such deviations are preceded by uniform conditions, it is seen 
that in all cases the other individual has maintained a single slope for high 
numbers of interactions (from 30-100 units, the number varying with 
each individual). If it be remembered that the duration of a slope is 
determined (within limits) by the ratio of its duration to the total duration 
(sum of both A’s and B’s slopes during that time), then it is evident that 
deviations from the fast-slow pattern follow after a slope has been main- 
tained longer than would have been the case if the rates of the two indi- 
viduals had adapted to one another. This adaptation is, of course, de- 
termined by the value of the function, log slope = k(ta/ta+ts)+constant, 
for the constituent individuals. 

Summary.—In recording the interaction of pairs of individuals, the 
measurements secured represent the durations of the separate actions of 
an individual alternated with his silences or inactions. If each five actions 
of an individual are summed and plotted cumulatively, a series of slopes 
is obtained, during each of which the individual maintains an approximately 
constant rate. 

The analysis of these slopes indicates that the duration of each slope 
is determined by the rate of the slope and its duration in chronological 
time (the sum of the durations of both A and B), the relationship having 
the form, log slope = k(ta/ta+ts)+ constant, within the limits of a band 
of constant width. It is suggested that the constant is invariant for the 
individual. 

Each individual is regarded as alternating between fast and slow rates. 
Although the values of these rates vary in a single conversation, the ratio 
of each fast rate to its accompanying slow rate is constant for the con- 
versation, but varies from event to event. 

Deviations from this pattern occur in more than half the cases. It is 
shown, however, that both representatives of the pattern (calculated) and 
its deviations vary in a uniform way. If the mean of the slopes of each 


> 
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conversation is calculated and its standard deviation, it is found that log 
a slope = k(log Mean slope) + constant, with all the points falling within 
a band of constant width. This relationship was the same for all observa- 
tions on four individuals. It was then assumed that deviations could not 
be due to individual differences in variability as measured by this function. 
It was then shown that deviations were in all cases preceded by slopes 
maintained over a period of time due to the operation of the suggested 
uniformity governing the duration of slopes at different values. 


1 Thomas, D. S., Loomis, A. M., Arrington, R. E., Observational Studies of Social 
Behavior, New Haven (1933). 

2 Chapple, E. D., Measuring Human Relations (to be published). 

3 Crozier, W. J., and Holway, A. H., Proc. Nat. Acad. Sci., 24, 3 (1938). 

4 Crozier, W. J., Determinisme et Variabilite, Paris, Hermann, 56 pp. (1935) (and 
citations to earlier work); Jour. Gen. Physiol., 19, 503 (1935-1936); Crozier, W. J., 
and Holway, A. H., Proc. Nat. Acad. Sci., 24, 3 (1938). 


ASTROPHYSICAL CONSEQUENCES OF METASTABLE LEVELS 
IN HYDROGEN AND HELIUM 


By O. Struve, K. L. G. HENYEY 
YERKES OBSERVATORY, UNIVERSITY OF CHICAGO 
Communicated January 9, 1939 


1. The metastability of an atomic level can influence its population only 
under conditions which differ materially from thermodynamic equilibrium. 
In or near equilibrium the populations of the levels are independent of 
transition probabilities and depend only upon the temperature of the 
system. On the other hand, when the deviations from equilibrium are 
marked, as, for example, in the case of a gas excited by diluted black-body 
radiation (nebulae, outer shells of stars), cyclical processes play a funda- 
mental réle. Observationally the metastability of a level produces dis- 
cernible effects when the interval of time between successive excitations is 
less than the lifetime of the level. Under these conditions, the metastable 
level can have a population of the same order of magnitude as in the case of 
equilibrium. Consequently, we may expect to observe in the continuous 
spectra of stars shining through the gas absorption lines arising from meta- 
stable levels, in addition to those arising from the ground level. Ofcourse, 
for sufficiently great thicknesses of the gas we may observe weak ab- 
sorption lines even when the metastability is not completely effective, 
provided that the number of atoms in the level is greater than about 10” 
in a column whose cross-section is 1 cm.? along the line of sight. 

The existence of metastable levels in hydrogen (2S) and helium (21S 
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and 2%S) is particularly interesting since the spectra of these elements can 
be observed astrophysically over an enormous range of conditions: from 
stellar atmospheres to-nebulae and, in the case of hydrogen, even to inter- 
stellar space. The emission lines of H and He are very strong in the Orion 
Nebula.! in absorption, only the He line \ 3889, which arises from the 
metastable level 2°S, is seen? superposed over the continuous spectra of 
stars which shine through the nebula. No H lines arising from the meta- 
stable 2S level have been found, in spite of the fact that H is almost cer- 
tainly much more abundant than He. It is the purpose of this paper to 
show that this apparent discrepancy may be explained by a change in the 
lifetime of the metastable H level which depends upon the electron density 
of the gas. 

For an excited level which connects directly with the ground level, the 
population in the presence of diluted radiation is: 


» (la) 


where W is the dilution factor. For a metastable level whose forbidden 
transition probability is A» the population is! 


hve 


(1b) 
An 

2. The Balmer lines are very strong in emission in the Orion Nebula, 
and intensity measurements by Ambarzumian! indicate that the numbers 
of atoms in several quantum levels higher than the second are all of the 
order of 104cm.~*. The number of atoms in the second level is not known, 
because we cannot observe the Lyman emission lines and because observa- 
tions fail to show any definite trace of the Balmer absorption lines. The 
latter are difficult to observe because all stars located in the Orion Nebula 
contain ordinary stellar absorption lines of hydrogen. However, in some 
of these stars the stellar Balmer lines are greatly broadened by the ionic 
Stark effect and by axialrotation. The absence of strong and sharp Balmer 
absorption lines of nebular origin is therefore well established. It is 
probable that if the number of atoms in the second level m,. = 10!*, the 
absorption line of the nebula would have been observed. We shall assume 
that m. < 10'*. In the case of hydrogen we may assume that vy ~ 143. 
Hence from (la) and (1b) 


A 1012 
(2) 


Accordingly 
Tes < 1 sec. (Orion Nebula) (3) 


We 
ng = 
= 
21 
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3. We have good evidence that interstellar space contains a consider- 
able amount of atomic hydrogen, and a rough determination of the 
intensities of interstellar emission lines in certain regions of the Milky 
Way gives mj = 10 cm.~*. This value is probably not representative of 
all regions of interstellar space, since it depends upon observations of 
selected regions in which the brightness of the emission lines is greater 
than the average. A better estimate for average conditions is 


= 1 


This is the number of third-level atoms in a column having a cross-section 
of 1 cm.” and extending to the effective limit of vision, which we estimate 
to be of the order of 1000 parsecs = 3 X 1074 cm. There is very good 
evidence from the spectrum of Nova Lacertae, whose distance is about 900 
parsecs,‘ and from other objects, that there are no observable absorption 
lines of hydrogen. Hence m; < 10! and by formula (2): 
“ < 10% and 725 < 10‘ sec. (Interstellar Space). (4) 
4. Bethe’ has discussed theoretically the question of the lifetime of 
the 2°S.,, level of hydrogen. He states that the lifetime due to spontaneous 
transitions to the ground level is of the order of several months. If, 
however, the atom is in an electric field the atom can also pass to the 
ground level by first going to the 2?P,,, level. Indeed, under most astro- 
physical conditions the fields, due to the presence of ions, govern wholly 
the lifetime of this level. Bethe shows that, for a weak field, the ratio of 
the lifetime of the 2.S level to the lifetime of the 2P level is given by 


(22) (5) 


T2P k 


where 1/48 is the lifetime of the 2P level and k/2m is the Stark splitting 
of the 2S level, that is, 


Here F is the electric field expressed in volts per centimeter and x25 is 
the matrix element of a coérdinate of the electron in the hydrogen atom. 
The average field due to the presence of electrons has also been con- 
sidered by Bethe. If ” is the number of electrons per cm.*, the electron 
occupies a volume 1/n which, if replaced by a sphere, has a radius 


(7) 
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Bethe takes one-half of this radius as the mean distance between the 
hydrogen atom and the nearest electron. The field of the latter is, there- 
fore, 


9\2 
F = 300e (?) volts/cm. (8) 


The fields due to other electrons can be neglected since they can be regarded 
as forming a uniform charge distribution. 

Equations (5) to (8) can be combined so as to express the ratio (5) in 
terms of and of atomic constants. Inserting the numerical values of 
the constants,® we find: 

= 2.3 X (9) 


T2P 


1 
Since? — = 48 = 6.3 X 107-3, we have 
ms = 4 X 10'n-*’?, (10) 


It is also of interest to investigate the order of magnitude of this effect 
for the metastable level 24S of helium. In a field the 21'S level interacts 
with the 2'P level and becomes partly a P level. The ratio of the life- 
times of the 2'P and 21S levels are in the ratio of the Stark splitting to the 
separation of the interacting levels.* Using hydrogenic wave functions, 
Foster® has computed the Stark splitting for a weak field as 0.0285 F? 
cm.~', where F is expressed in 10° volts/em. The separation of the levels 
is 5857 cm.~'. Hence, expressing F again in terms of ”, we have 


Tas —4/3 


A comparison of this equation with (10) shows that the effect under con- 
sideration is of a completely different order of magnitude for the two 
cases. The importance of the effect for hydrogen is due to the fact that 
the 2S and 2P levels have exactly the same energy. For the metastable 
2%S level of helium the effect of a field is even smaller than for the 2'S level. 
Actually, under most conditions, the lifetime of the 24S level will be gov- 
erned by other factors, particularly by collision transitions to the 23S level. 

5. For interstellar space we have, approximately :* 


Ne = 30 cm.~*. 
Hence, formula (10) gives 
72s = 4 X 10° sec. (Interstellar Space). (11) 


Before we compare this value with the observational result (4), we must 
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multiply the latter by 4 to take account of the fact that the levels 2P,/2 
and 2P3 . are not metastable. We conclude that the theoretical lifetime 
is consistent with the observational value. 

For the Orion Nebula n,; is not known. Hence, we determine », from 
(10) in such a manner that (3) is obeyed: 


Ne > 2X 108 cm. 


We have no means to verify this result by direct observation. For the 
Trifid Nebula, in which the density is probably considerably smaller than 
in the Orion Nebula, Rosseland” found 


N= 


The somewhat higher value for the Orion Nebula appears entirely plausible. 

6. O. C. Wilson? has found that in the Orion Nebula He I (2%S-3'P) 
d 3889, appears as a strong absorption line. No accurate determination 
of the equivalent width is available, but we estimate m. = 10'* cm.~?. 
The emission lines of He I are strong, but no determinations of m3; have 
been made. Perhaps we may estimate n; = 10*cm.~*. This would give 
from (2): 


— = 10” and ta5 = 10? sec. (Orion Nebula). (12) 


Accordingly, the metastability of the He I 2*S term in the Orion Nebula 
is much more pronounced than that of the H 2S term. This is in good 
agreement with the theory because we have seen that the theory would 
give an exceedingly long lifetime for the term. 

The question arises whether the astronomical evidence is compatible 
with a practically infinite lifetime of the 2°S level, as has been suggested!! 
from physical theory. Apparently this is not the case. Formula (1b) 


A A 
holds if W < , we have approximately: 


31 31 


hv 
= me *T. (13) 


If this were true (13) and (la) would give 


W = 
It is probable, however, that for those parts of the Orion Nebula to which 
the observations refer W is more nearly of the order of 10~'* or 10~-". 
Hence, we should be inclined to regard (12) as real. 
The He I line (2'S—4'P) d 3965 has been looked for in the spectra of the 
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same Orion stars which show strong (2%S-3*P). No absorption line has 
been found. However, since all lines of the singlet system of He are 
much weaker than the corresponding lines of the triplet system, we cannot 
assert, from the observational evidence, that the lifetime of the 24S level 
is necessarily shorter than that of the 2°S level. A careful search for the 
line \ 5016 should settle this question. 

7. The spectrum of Nova Lacertae shows conclusively that for inter- 
stellar space the He I line (2°S-3*P) gives m. < 10'*. Since no emission 
lines have been observed we cannot determine the upper limit of As/Az. 
The ionization of the interstellar gas is much lower® than that of the Orion 
Nebula; hence the absence of the He I lines is not surprising. 

8. Formula (10) shows that for relatively low electron pressures the 
2S state of H completely loses its metastability. Thus, when 


10" cm.-* 


the lifetime becomes of the order of 10-* sec. This is of the order of the 
electron density in the chromosphere. In outer shells of B and O stars, 
where the electrons come almost entirely from the ionization of hydrogen, 
we estimate m, = 10° cm.~*. Hence the hydrogen lines should behave 
as though they originated from a metastable level.'!2 However, it is 
possible that in some relatively dense shells n, becomes large enough to 
reduce the lifetime of the 2S state of H to such an extent that the Balmer 
lines will lose with respect to lines of other elements which originate from 
the ground level or from a metastable level. One is reminded of the great 
strength of the He I line \ 3889 in y Cas'* which surpasses the Balmer 
lines, although in normal giant spectra (where Stark effect is absent) the 
Balmer lines are always much stronger than the He lines. High density 
and strong ionization reduce the lifetime. Hence, in late-type giants, 
conditions are favorable for a high degree of metastability of the 2S level. 
The enhancement in @ Ori of the multiplets (a5D-2’D°) and (a5D-z’P°), 
which originate from the ground level of Fe I, show that departures from 
thermodynamic equilibrium are conspicuous. It is tempting to explain 
the great strength of the Balmer absorption lines in supergiants of late 
spectral type as a consequence of the metastability of the 2S level. How- 
ever, it is premature to discuss this matter further until accurate measure- 
ments of the Paschen lines have been made. 

1V. Ambarzumian, Zeits. Astrophys., 6, 112 (1933); Poulkovo Obs. Circ., 6 (1933). 
This expression is valid only when W is small compared to Ao /Au. 

20. C. Wilson, Pub. Astr. Soc. Pacific, 49, 338 (1937). 

3Struve, Proc. Nat. Acad. Sct., 25, 36-43 (1939). 

4J. A. Pearce, Pub. Astr. Soc. Pacific, 49, 149 (1937). 

5 Handbuch der Physik, 24, Pt. 1, 452 (1933). 

6 The matrix element x**:p is 1/3 times the Hartree unit of length, as may be verified 
from the formulae 34.2, et seq., of Bethe, into which we insert m = 2,7 = !/,andm = 1/2. 
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7 Bethe, op. cit. Table 17, page 144. 

8 Condon and Shortley, The Theory of Atomic Spectra, p. 413 (1935). 

9 Proc. Roy. Soc., A117, 137 (1937). 

10 Theoretical Astrophysics, p. 319 (1936). In his earlier book, Astrophysik auf Atom- 
theoritischer Grundlage, p. 235 (1931), Rosseland gave n, = 470 cm.~%. 

11 We are indebted for stimulating discussions of this point to Dr. D. H. Menzel of 
Harvard University and Dr. G. Breit of the University of Wisconsin. 

12 Struve and Wurm, Astrophys. Jour., 88, 107 (1938). 

13 R. Baldwin, Astrophys. Jour. (in press). See spectrum (a) of his illustration. 


THE RESOLUTION OF SIX TESTS INTO THREE GENERAL 
FACTORS 


By EpwIn B. WILSON AND JANE WORCESTER 
HARVARD SCHOOL OF PUBLIC HEALTH 
Communicated January 14, 1939 


It has been pointed out that if ” tests are to be resolved into / generals 
and n specifics the expression 


(n+) 


must be zero or negative and that when T = 0 the resolution if possible 
at all is unique in the sense that no arbitrary parameter is involved in the 
solution but not necessarily in the sense that there shall be only one solu- 
tion.! The case of six tests and three factors is the simplest in which the 
possibility of multiple distinct solutions arises. It is the object of this 
note to exhibit a 6 by 6 correlation matrix which has two solutions for the 
resolution.? 
The problem is to determine six communalities h? so that the matrix 


hi "12 113 "14 "15 "16 
113 hs 134 135 36 
T1400 hi 146 


shall reduce to rank 3. There is no four-rowed minor which does not 
contain at least two communalities; consider therefore 


N14 "15 "16 

125 126 0 
2 

h3 134 135 136 

hk 


73 
. 
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This expands into a quadratic equation in the communalities h} and h{ as 


hh? N16 "14 N16 Nis Ne "13 Nia Nis 16 
34 2 
O 0 135 136 O 135 36 


O 


Including the above, three equations may be obtained by non-equivalent 
permutations among the subscripts 1, 2, 5, 6; these three are linearly 
related, indeed if their signs be taken aright their sum is 0; of the set of 
three, however, any pair are usually independent and may be solved for 
h? and hj to give a pair of solutions, ordinarily different. When either of 
these pairs is inserted in the place of h} and hj, four-rowed minors may 
be found which contain only one other communality and hence the com- 
plete solution contains just two sets of communalities. 
The following matrix 


h? 56 .16 .48 .24 .64 
16 .20 kh 
4 18 & 


may be reduced to rank 3 by either of the following set of communalities: 


hy = .64 or h; = 1015 .425616 
1173 

5 2 = = 

he = .85 or hy 1300 . 902308 
311 

= .06 or = 4900 .063569 

= .56 hi = = .546923 

= .56 or = = - 2: 

hi = .50 or h§ = a 386667 

s= -50 or lis = 300 = - 7 
998 

.93 or = 1000 298000 


It will be observed that although some of the communalities in the two 
sets are nearly equal, others, notably hj and hj (and, presumably, because ' 
of their nearness to 1, h{) are decidedly different. It may further be noted 


| 
| 
| 
| 
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that in both sets some of the communalities are tolerably far removed from 
the value of the largest correlation coefficient in their row or column which 
is the provisional value assigned by Thurstone in his centroid method of 
determining the factor loadings.* 

One set of factor loadings for the two solutions is: 


I II III II’ III’ 
x 8 0 0 6524 0 0 
7 0 .4068 0 
X3 2 1 2453 — .0259 .0515 
6 ~ 2 7358 .0265 
X5 3 4 — .1530 
8 5 .9810 . 1835 


It will be noted that the factor loadings are different, often decidedly 
different. Furthermore, although the set I, II, III or the set I’, II’, 
III’ may be changed by an orthogonal transformation of * degrees of 
freedom operating upon them, the two sets cannot be obtained one from 
the other by any orthogonal transformation operating upon them*—the 
solutions are distinct. 

If we should take the correlation matrix 


h? 


and try to determine the communalities so as to reduce the rank of the 
matrix to 3 we should find only one set of solutions all of which lay between 
0 and 1 and hence were possible; the other set, though real numbers, would 
contain at least one communality which was either negative or greater 
than 1. The solution for the possible set of communalities and of factor 
loadings is 


he 644 0 
0 
&.-2 
he 4 
he “a. 4. <2 


= 
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If we should take the correlation matrix 


he 56 .16 .48 .24 .64 
.20 .66 .51 .86 
16 .20 kj 18 .07  .25 
48 66 .18 hi .76 
8 2 .% 


we should find that no set of communalities could be obtained which would 
reduce the rank to 3 and hence permit a resolution into 3 generals and 6 
specifics. The three matrices do not look notably different. 

If the pentad expression’ be defined: as 


— 14103735045 — 113% 15% 237%24%45 — 1127157 2473435 
— 1314124705735 — — 1147157 23725734 


the quadratic equation for any communality, say 3, becomes 
+ + = 0, 


where H and K, like P, are skew symmetric in the letters 1, 2, 4, 5, 6 and 
are in fact as follows: 


2 2 2 2 2 
K3.10456 = 30456 + 134235612 + 36124 + 
= 151 [T4556 + 1467 45( 116725 = 115726) + 


+ nine similar terms in which the pair 1, 2 1s replaced by any other 
pair of two from the five numbers 1, 2, 4, 5, 6. 


The analytical conditions for resolubility into two sets or just one set or 
no set of three generals and six specifics are naturally complicated; it 
is probably better to proceed to effect a solution arithmetically. 

If one should go up to the next case in which the resolution may be 
“unique,’ viz.,n = 10,/ = 6, there may well be a higher multiplicity than 
two for distinct solutions. 

If one should undertake to find solutions by Thurstone’s centroid method 
it is clear that some method of successive approximations and convergence 
would be necessary—and necessarily very tedious. It may well be that 
the actual determination of any or all solutions in a particular case would 
be of no value, but it does seem to be.of some value to recognize that the 
method of factor analysis as a mathematical method admits of a plurality 
of solutions in some cases. In so far as three general factors are deter- 


VoL. 25, 1939 PSYCHOLOGY: WILSON AND WORCESTER 77 


minable from six tests there may logically be two distinct solutions—and 
therewith two distinct sets of specifics but the team tests for corresponding 
specifics in the two sets must be the same.’ 


1E. B. Wilson, Jour. Gen. Psychol., 2, 153-172 (1929), especially p. 156. 

2 The authors are indebted to the Carnegie Corporation for a grant to cover the cost 
of the calculations involved. 

3L. L. Thurstone, The Vectors of Mind, p. 89. As a matter of fact hj in either set is 
less than any correlation coefficient in its row or column, and h? is greater than any. 

4 Of course, as the two sets are each three mutually orthogonal unit vectors, either 
can be carried into the other by some orthogonal transformation operating not in the 
3-space defined by either set but in the 6-dimensional space defined by the six vectors 
of the two sets—this being not the 6-space defined by the 6 tests. 

5 T. L. Kelley, Crossroads in the Mind of Man, p. 58. 

* If one should apply Thurstone’s method without proceeding to successive approxi- 
mations one would find a set of generals, with their factor loadings, which was not 
identical with either of the two sets I, II, III or I’, II’, III’. If it were true that the 
set found could be rotated into either of these sets by an orthogonal transformation 
operating upon it, the lack of identity might be of no significance by virtue of Thur- 
stone’s assumption that in any case the set found is without psychological significance 
until it has been so rotated as to be psychologically interpretable; but as the solution 
found by Thurstone’s method is not in general in the 3-space of I, II, III or of I’, II’, 
III’ no rotation of it in the sense in which he uses that word will bring it into agreement 
with either of those sets which the general theory indicates might be of significance. It 
may well be that the Thurstone centroid solution and the true distinct solutions I, II, 
Ili and I’, If’, III’ all lie nearly enough in the same 3-space to be practically indis- 
tinguishable when one has regard to the sampling errors of the correlation coefficients 
determined from relatively small samples. It may also well be that within the freedom 
permitted by those sampling errors the question as to whether the resolution can be 
made in two, one or no ways is itself indeterminate. If more than 6 tests are to be 
resoluble into three general factors plus specifics very complicated analytical conditions 
must be satisfied by the correlation coefficients. If 6 tests are resoluble into three 
general factors plus specifics in two ways one may conceivably build up additional tests 
which would be consistent either with the solution I, IJ, III or with the solution I’, 
II’, III’, but ordinarily not with both. We may thus have two sets of m tests, m > 6, 
resoluble into three generals and m specifics, the generals as well as the specifics being 
different, even though there be six tests common to the two sets. 

7E. B. Wilson, Proc. Nut. Acad. Sct., 20, 193-196 (1934). 
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TEMPERATURE AND THE CRITICAL INTENSITY FOR 
RESPONSE TO VISUAL FLICKER. II 


By W. J. CRozIER 
BIOLOGICAL LABORATORIES, HARVARD UNIVERSITY 


Communicated January 13, 1939 


1. Initial tests of the relation between temperature and the curve of 
flash-frequency (F) vs. flash-intensity (J) critical for response to flicker! 
showed that change of temperature merely shifts the F — log J contour 
on the log J axis, without essentially altering its form. Rise of 
temperature moves the curve to lower critical intensities. This fact is 
of primary significance for the theory of the flicker-response contour.” 
The exact mode of dependence of the shift in position upon temperature, 
as revealed by more elaborate experiments with various animals, thus 
assumes particular importance. The first experiments! were made at 
3 temperatures only. For an insect (Anax, nymph) and a teleost (sun- 
fish, Enneacanthus) it appeared that a somewhat unusual relationship 
obtains. Theory calls for the reciprocal of the critical intensity at a fixed 
flash-frequency as the proper measure of excitability. The temperature 
characteristic® (u in the Arrhenius equation) for 1/J was not constant but 
increased with rise of temperature. Hence it was suggested! that visual 
excitability measured in this way might well be so complexly determined 
as to make its analysis quite involved. 

More complete investigation was made of the flicker-response contour 
of the turtle Pseudemys.* Here a simplex curve is found, without compli- 
cations due to the gross morphology of the eye (as in most arthropods*) and 
without the involvement of two sets of visual response-elements (‘‘rods”’ 
and ‘“‘cones’” as in most vertebrates). As fundamentally in all other 
cases investigated, the F — log J contour is a probability integral. The 
adequacy of this formulation arises directly from the fact that a large 
number of elements are available for the determination of the response, 
and that each of these fluctuates in its contribution. Elevation of tem- 
perature reduces the critical intensity for each of these elements in the 
same proportionate way, if the population is homogeneous with respect 
to type of chemical organization. Hence if excitability is governed in 
the fashion typical for physiological processes involving speeds or fre- 
quencies,® we expect 1/J (F constant) to obey the Arrhenius equation, 
with constant w’s. For Pseudemys this was found.‘ 

The exceptional instances first encountered were therefore reinvestigated 
with especial care. Processes demonstrably complex in their quantitative 
dependence on temperature should provide significant tests of the general 
theory of temperature characteristics. Of these the most important for 
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the view that the control of speeds and frequencies in biological systems 
is essentially simple rather than complex, are those instances in which u 


increases steadily with rising temperature. 


Processes of this kind are 


necessarily governed by:the resultant of 2 or more different concurrent 
reactions with different ‘‘activation energies” simultaneously contributing 


to the control of the occur- 
rence of the observed result.’.* 
Such processes have seldom 
been established in biological 
observations—a fact which is 
in itself significant. The mean- 
ing of the data in the single 
most fully illustrative case, 
involving the respiration of 
sea urchin eggs,’ has been 
completely missed.!° 

2. Determinations of criti- 
cal intensity for response of 
Anax nymphs to visual flicker 
have been made at 2 flash- 
frequencies (F = 20/sec. and 
F = 55) at temperatures be- 
tween 8° and 36°.!1 The 
technique, methods of calcu- 
lation and tests for statistical 
homogeneity of the material 
have been described  else- 
The form and 
maximum of the F — log J 
curve is independent of tem- 
perature. (A minor compli- 
cation not significant for the 
present discussion occurs at 
36° for F = 55 only.) The 
apparent temperature charac- 
teristic for 1/J is the same 
for F = 20 and F = 55, and 


a. 
10/ Tas 


FIGURE 1 

1/I,, at F = 20 (solid dots) and at F = 55, at 
various temperatures, for Anax nymphs, with 
light-time fraction = 0.50. Each JI,, is the 
mean of 30 measurements.'! Values for F = 
55 are multiplied by antilog 1.205. The excep- 
tional points at F = 55, 35°8 are by experi- 
mental analysis known to be due to special 
effects at this level;!! J,, departs by only 0.03 to 
0.09 log unit. The upward concavity on the 
Arrhenius grid shows that u in the equation 1/I 
= exp. (—u/RT) increases continuously as T 
rises. Consequently at least 2 different proc- 
esses participate concurrently in the control of 
1/J. The curvature found is accurately re- 
produced by the summation of the two proc- 
esses whose velocity curves are shown below, 
one with » = 19,200 and the other with »p = 
3400; their velocities are equal at 15.9°. 


increases steadily as ¢° increases (Fig. 1). 

The upwardly concave curve in figure 1, with log (1/Z) plotted against 
1/T°4ss., is easily accounted for as the resultant of the concurrent action 
of 2 processes with » = 19,200 and » = 3400, respectively, with velocities 
equal at 15.9° (Fig. 1). The suggestive correspondences of » = 19,200 
for dehydrogenations!! (as of sugars) and such processes as diffusion for 
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the low u, point the way for specific attempts at the identification of these 
processes through suitably designed experiments. 

3. The flicker contour for the sunfish Enneacanthus is found to present 
a different kind of situation.'* Its elucidation again illustrates the neces- 
sity for closely spaced observations of known relative precision for the 
interpretation of temperature functions. With Enneacanthus there 
occurs a break in the 1/J curve for a given F, at ca. 20°C. Above this 
temperature u is higher than from 8° to 20°. Relations of this type are 
uncommon, but not unknown,*'* and there is nothing in theory which 
forbids them. 
| One flash-frequency (F = 4) 
chosen in the so-called 

“rod” section and another 

(F = 25) in the ‘‘cone” part of 
the sunfish flicker contour." 
| The shape of the F — log I 
| curve is independent of tem- 
perature, in each of these seg- 
ments. But the break in the 
temperature curve (Fig. 2) is 
unmistakably emphasized by 
the fact that for the “rod” 
- part of the function the curve 
20. 50 is definitely shifted downward 
on the 1/J axis to a lesser 


7 . extent than for the cone part, 

1/I,, for F = 4 and F = 20, the former di- hal 20° Thi k it 
vided by 23,950, with the sunfish Enneacanthus. 
Each point is the mean of 30 observations. obvious that the critical tem- 


Points with tags are from an early series! For perature at ca. 20° is a real 


the “rod” (F = 4) and the “cone” (F = 20) phenomenon in this system, 
parts of the F = log J curve the y’s are identical. even if it were not already 


The critical temperature (20°) is marked by a 5 ; 
wider dislocation of the F = 4 measurements. clear that no single contint- 
ous curve will adequately 


describe the data at F = 25 (Fig. 2). 

The values of u» for the sunfish data are 8100, 14,400. The same values 
hold for both ‘‘rod”’ and ‘‘cone”’ sections of the curve. At any temperature 
the values of 1/J,, for F = 4 and F = 25 are in the ratio of 24,000: 1 or 
(below 20°) more. This is a rather striking illustration of the fact that 
uw for a biological system is essentially independent*** of the velocity 
of the underlying process (here measured by 1/J). 

4. That diverse magnitudes of » are encountered for flicker-excitability 
in different animals naturally signifies that the specific chemical control 
of the excitability is different. This does not necessarily mean that the 
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substances undergoing chemical change as result of illumination are 
different. Substrate materials in the excitability reactions could be the 
same, but catalysts different. 

It is just as likely, however, that the substances involved could even 
be the same, and their metabolic interrelations, but that in different kinds 
of animals different steps or links in the system are found in the slow 
“pacemaker” réle.4 The data on experimental modification of yu!” make 
this entirely possible. Differences in the locations of critical tempera- 
tures (Pseudemys, 20°; Enneacanthus 30°) prove that such differences in 
organization are real. 

Application of the hypothesis'* that the data of visual flicker are to be 
interpreted in terms of retinal photochemistry has led’* to the conclusion 
that for Enneacanthus there must be deduced a different photochemical 
stationary state equation for the ‘‘rod”’ and “cone” parts of the contour. 
For the “rods,” light and dark processes are each of second order, for the 
“cones,” light process is first order, dark process second. The equation 
does not in fact describe the F — log J curve,”® and is moreover incon- 
sistent with the direction of the temperature shift; apart from this, the 
invariance of » for “rods” and ‘‘cones’”’ (Fig. 2) shows directly that the 
chemical control of excitability is the same in both. 

This is not in conflict with the existence of 2 populations*! of sensory 
effects, as proved by their behavior in inheritance”? and by the present 
results with Enneacanthus. Since there is every reason to regard these 
2 populations of elementary sensory effects as resident in the central 
nervous system,”* there is a general basis for expecting their chemical 
organization to be of the same kind. That it is of the same kind, in 
Enneacanthus, is shown by the fact that with increasing proportion of 
light-time in the flash-cycle the rates of increase of the abscissa of inflection 
in the “‘rod’’ and ‘“‘cone”’ parts of the F — log J curve are identical.” 
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EXPRESSIONS FOR THE CURRENT IN THE BLOCH APPROXI- 
MATION OF “TIGHT BINDING” FOR METALLIC ELECTRONS 


By J. BARDEEN AND J. H. VAN VLECK 
UNIVERSITY OF MINNESOTA AND HARVARD UNIVERSITY 
Communicated January 14, 1939 


Various approximate methods have been used for the calculation of 
metallic wave functions.' Of these, the most important are: (1) the 
approximation of “tight binding” in which it is assumed that the atoms 
are relatively far away from one another, so that the overlap of the wave 
functions of the electrons in neighboring atoms is small;? (2) the opposite 
limiting case, the approximation of nearly free electrons, in which it is 
assumed that the variations in the potential energy are small in comparison 
with the kinetic energy of the metallic electrons’ and (3) the method of 
Wigner and Seitz‘ and its extensions by Slater,® which is the most accurate 
of all. We wish to discuss the expressions for the current in the first of 
these methods. In this method, the Bloch wave functions are approxi- 
mated by a linear combination of atomic wave functions. 

The Bloch form for the wave function of a metallic electron is y, = 
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exp(7k-r)u, (r), where u,(r) has the translational period of the crystal 
lattice. The vector k is the propagation vector of the electron wave, and 
is in magnitude 27/\, where ) is the wave-length. The standard expression 
for, say, the x-component of the current carried by an electron in this 
state is: 


i, = = —(eh/Ampt) S Ox) — Ox) (1) 
The Bloch expression for the same component of the current is® 
i, = —(2ne/h)(OE/0k,) (2) 


where £ is the energy, expressed as a function of the propagation vector 
k. These two expressions are, of course, equal if Y, is an exact solution of 
the Schrédinger equation: 

— Ay, + = (3) 
in which V(r) is the periodic potential of the crystal lattice. If yy is 
determined by some approximate method, there is no reason to suppose 
that (1) and (2) will give the same result. In fact, if hydrogen-like atomic 
wave functions, e~“, are used for the calculation of the Bloch wave func- 
tions in the approximation of ‘‘tight-binding”’ it is found, as will be shown 
more explicitly later, that (1) is just one-third of (2).’ 

Such disagreement on the current is perhaps at first surprising, because 
one might expect that, at least in the limiting case when the distance 
between the atoms is very large, the approximate wave function so derived 
will be very close to the exact wave function, obtained by solution of the 
Schrédinger equation. However, in this limit the current goes to zero, 
and while both (1) and (2) approach zero, the ratio between them does 
not approach unity. 

It is the purpose of the present paper to point out this discrepancy, 
explain somewhat its causes and to examine the type of wave function 
for which exact agreement is to be expected. 

Let us first examine the nature of the field in which the electron moves. 
We may assume that the potential in the crystal lattice is fairly constant 
in the regions between the atoms, and, without loss of generality, this 
constant potential may be taken equal to zero. If the potential in the 
neighborhood of any atom is U(r), the total potential in the crystal is* 


V(r) = 2, U (|r — (4) 


where the sum is over all lattice points (designated by the vectors p;). 
Here U(r) cannot be taken to be the true atomic potential (i.e., the field 
of the ion about which the valence electron of the free atom moves), as 
is generally assumed, because the latter will fall off as 1/r for large r, and 
the sum (4) will be infinite for an infinite lattice. Rather, one should use 
a screened field which in the neighborhood of any atom is very nearly equal 
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to the field of the ion, but which falls off exponentially instead of as 1/r 
at large distances. The screening of course is due to all the valence 
electrons except the one under consideration. These considerations 
indicate that to compute the Bloch functions one should use, instead of 
the true atomic wave functions, functions which are determined from a 
screened ionic field. We show below that if this is done the expressions 
(1) and (2) give the same value for the current, at least in the limiting 
case when the atoms are far apart. The disagreement on the current 
mentioned above, which occurs if true atomic wave functions are used, 
thus seems to be connected with the choice of the crystalline potential. 

We proceed with the mathematical proof of the foregoing statements. 
The equation for the ‘‘atomic’’ wave functions ¢(r) from which we will 
construct the Bloch wave function is 


— Ag(r) + U(r)g(r) = Eng(r). (5) 
For simplicity we assume that ¢(r) is a radial (s) function. We also 


assumed that U(r) falls off much more rapidly than 1/7 for large r, so that 
the asymptotic form for ¢(r) which holds if U(r) = 0 is: 


g(r) ~e (6) 
in which a is determined from 
Ey = (7) 
The true atomic wave function would instead have the asymptotic form 
e “, which is appropriate for an electron moving in a coulomb field. 


The wave function of an electron moving in the crystal lattice is formed 
from an appropriate linear combination of the ‘‘atomic’’ wave functions.® 


ve = exp[ik-p;] o(\r — ol), (8) 
the sum running over all lattice points. 
In this approximation the energy, E,, is:* 
By = Ey + — Vie) — UG) e(*)dr. (9) 
The Bloch expression for the velocity of an electron in the state k is 
= (2n/h)(OE/Ok,) = —(2x/h)z; exp[—ak-p,], (10) 
in which J, is the integral: 
Selle — Ve) — ur) 
= — — pml) )o(r)ar. (11) 
m + 0 


The second form follows from the first by substitution of the expression 
(4) for V(r). Since we are assuming that U(\r — p,,|) is small except in 
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the neighborhood of p,,, the only significant term in the sum over m is 
that for which m = 1, so that 


— pl) U(r — m|)e(r)dr. (12) 
Using the Schrédinger equation (5) for g(r), we see that J; may be written 
in the form: 
1, = + — (13) 
Since 
Ag(|r — pil) = — ail), (14) 


(the second Laplacian is taken with p, as the variable) we may express 
I, in terms of the overlap integral 


g(o) = So(lr — pl)y(r)dr. (15) 
Thus 


I; = A, + Eo} 


h? 2dg 
-{ + - 


dp? + Eat} 


pdp (16) 


e=p, 


To obtain the second torm we have made use of the fact: that 9(r) and 
consequently g(r) are radial functions. 

The ordinary expression for the velocity is obtained by substituting 
the expression (8) for y,, into equation (1). We then have: 


v, = (th/2ru)z, So(r)(Oe(\r — py|)/Ox)dr. (17) 
Now 


— p;|)/Ox = —Oe(|r — p|)/O(m)., (18) 
so that 


— 
— (01) (19) 


In view of (16), the two expressions (10) and (19) for the velocity will be 


equal only if 
+ pdp \p dp/’ (20) 


or if (cf. Eq. (7)) 


Ux 


_ — = 0. (21) 


= 
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Thus we must have 
g=c” (22) 


where c is a constant. 

The overlap integral g(p) will have this form if the major contribution 
to the integral comes from the region where the potential is fairly constant 
(and, by assumption, equal to zero) because in this region the wave function 
will have the form ofr. It is not difficult to show that the overlap 


ealt—al 
integral, —— dr, is proportional to and thus has the 


desired form. 
On the other hand, for hydrogen-like wave functions of the form e~°’, 
the overlap integral is proportional to® 


e {1 + ap + 1/sap)*}, 


and with tnis expression (19) is just one-third of (10). 


1 See, for example, M. F. Mott and H. Jones, Theory of the Properties of Metals and 
Alloys (Oxford, 1936) Chap. II, or A. Sommerfeld and H. Bethe, Handbuch der Phystk, 
2nd Edition, 24, 2, p. 370, Berlin (1934). 

2 F. Bloch, Zeits. Physik, 52, 555 (1928). 

3R. Peierls, Ann. Physik, 4, 121 (1930). 

4E. Wigner and F. Seitz, Phys. Rev., 43, 804 (1933). 

5 J. C. Slater, Jbid., 45, 794 (1934). 

6 Cf. Mott and Jones (reference 1), p. 93. 

7 Analogous considerations apply to the corresponding two ways of computing the 
diamagnetic current induced by a magnetic field. [F. London, Jour. de Physique, 8, 
397 (1937).] If the diamagnetism is computed directly from an integral similar to (1) 
rather by differentiation with respect to the field strength [the analog of (10)] and 
if hydrogenic wave functions are used, London’s estimate of the resonance integral 
necessary to account for the diamagnetism of aromatic compounds is raised from 4.4 
to 13.2 electron volts. The fact that the latter value is unreasonably high need not 
cause concern, as the hydrogenic model is not really applicable. Anyway, the change 
by a factor of three cancels out in computing the relative amounts of diamagnetism of 
different aromatic substances, so that the ability of his theory to predict the variation 
of diamagnetism with chemical composition is not affected. 

8 Cf. Mott and Jones (reference 1), p. 65 and following. 

9 Cf. W. Heitler and F. London, Zeits. Physik, 44, 455 (1927). 
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THE OXIDATION-REDUCTION POTENTIALS OF LAKE WATERS 
AND THEIR ECOLOGICAL SIGNIFICANCE 


By G. E. Hurcuinson, E. S. DEEVEy, JR.,1 AND ANNE WOLLACK 
THE OSBORN ZOOLOGICAL LABORATORY OF YALE UNIVERSITY 


Communicated January 3, 1939 


Since the classical studies of Thienemann,? it has been realized that the 
nature of the profundal Chironomid fauna of a lake provides a valuable 
criterion in lake classification. It has usually been believed that oligo- 
trophic lakes with a rich supply of hypolimnetic oxygen have a Tanytarsus 
fauna, eutrophic lakes, with a poor supply, a Chironomus fauna. Lund- 
beck* elaborated this generalization, characterizing various mesotrophic 
categories, without adding much information as to the chemical condi- 
tions of the different subtypes. The work of Miyadi‘ indicates that 
inconsistencies, notably involving ‘‘oligotrophic Chironomus’’ lakes, may 
arise in some regions. In the course of a study of the benthic faunae of 
Connecticut lakes, one of us (E. S. D.) has recognized comparable in- 
consistencies, as for instance when Lake Quassapaug, Middlebury, is 
compared with Job’s Pond, Portland. At the end of summer stagnation 
the oxygen content of the bottom water of the former lake is, in spite of its 
greater depth, lower than that of the latter, yet Quassapaug is a typical 
Tanytarsus lake, while Job’s Pond resembles Miyadi’s ‘“‘oligotrophic 
Chironomus”’ lakes. Chemical studies (Hutchinson®) have indicated large 
amounts of ferrous iron in the hypolimnion of Linsley Pond, a typical 
eutrophic Chironomus lake, while Jahn® found that the bottom water of 
Quassapaug has a redox potential but little lower than that of the surface 
water; comparison of Yoshimura’s’ figures for hypolimnetic iron in some 
of the lakes studied by Miyadi strengthened the suspicion that oxidation- 
reduction potentials might be of interest in relation to the Chironomid 
typology of lakes. We have accordingly measured the potentials set up 
on a bright platinum electrode immersed in freshly collected lake waters, 
and at the same time have determined ferrous or total iron, total manga- 
nese, pH and oxygen. All potentials were determined between 15° and 
20°C.; to the nearest centivolt, they may be taken as correct for 18°C. 
Hill’s’ a-a’ dipyridyl method was used in the determination of ferrous 
iron. 

In general the freely circulating epilimnetic waters give potentials of 
from 0.46-0.50 volt, referred to the standard hydrogen electrode; such 
values are to be expected in well-oxygenated waters in the absence of 
reductants (Cooper’). Two lakes give somewhat lower surface values 
(Queechy, Canaan, N. Y., 0.43 volt, and Wononskopomuc, Salisbury, 
0.44 volt). In the hypolimnion of most lakes the potential is found to 


: 

: 
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vary with the amount of ferrous iron present, though owing to the high 
dilution and complex composition of the water theoretical treatment can- 
not be attempted. In only one case, that of the bottom water of Queechy, 
are both the redox potential and ferrous iron (Eh 0.15 volt, Fett 0.4 


mgms. per liter) so low as to suggest the presence of a non-ferrous poising 
system. 
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FIGURE 1 


Data for Lake Quassapaug (Q), 3 Oct. 1938 and for Highland Lake (H) 7 Oct. 
1938. 


In figure 1, data for Lake Quassapaug, with an oligotrophic Tanyiarsus 
fauna, and for Highland Lake, Winchester, Connecticut, with a mesotrophic 
Endochironomus fauna are set out, while in figure 2 data for the eutrophic 
and “‘oligotrophic’”’ Chironomus lakes, Linsley Pond, North Branford, Con- 
necticut, and Job’s Pond are presented. It will be seen that the Eh 
curve shows little variation in the 7anyiarsus lake, slightly greater hypo- 
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limnetic reduction in the Endochironomus lake, a strong microzonal fall 
in potential in the oligotrophic and marked reduction throughout the 
hypolimnion in the eutrophic, Chironomus lake. All have ‘‘eutrophic”’ 
oxygen curves. In certain mesotrophic lakes Chironomid larvae are 
found which, though yielding on emergence adults of a new species of 
Chironomus (s. str.) near hyperboreus, have but one pair of, or no, ventral 
blood gills, a condition common in mesotrophic species. This form is 
provisionally designated as ‘‘mesotrophic Chironomus.” It occurs with 
Tanytarsus in Lake Wononskopomuc, which shows a very slight fall in 
potential, to 0.41 volt, at 30.2 m., and without Tanylarsus in East Twin 
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FIGURE 2 
Data for Job’s Pond (J), 4 Oct. 1938 and for Linsley Pond (L), 16 Sept. 1938. 


Lake, Salisbury, where the potential falls microzonally to 0.33 volt at 
22.5 m. Mount Tom Pond, a true Chironomus lake intermediate in gen- 
eral character between Linsley and Job’s Pond, shows a fall from. 0.46 
volt at 8 m. to 0.26 volt at 13 m. Queechy Lake, a rather oligotrophic 
true Chironomus lake, shows a microzonal fall to 0.15 volt at 0.5 m. above 
the deepest mud, in this resembling Job’s Pond; alone of all the lakes 
studied a low potential occurred without much ferrous iron being present 
and a sulphide poising system may well be involved. It would appear 
therefore that where the water approximately 0.5 m. above the deepest 
mud has a potential of above 0.4 volt, the mud itself supports a pure or 
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mixed Janylarsus fauna, where below 0.3 volt a true Chironomus (four 
ventral blood gills) fauna is found, while in the intermediate range either 
Endochironomus or mesotrophic Chironomus, with a reduction of ventral 
blood gills, occurs. It is important to note that muddy water from 
Wononskopomuc shows that very low potentials (0.15 volt) are prevalent 
in the mud itself in the vicinity of the animals, so that the relationship here 
established is primarily of typological interest, indicating the reducing 
power of the mud as it affects the open water, rather than of significance 
in the respiratory physiology of the mud-inhabiting animals themselves. 

The general data incorporated in the above account were partly ob- 
tained from a study supported by the BAcHE Funp of the NATIONAL 
ACADEMY OF SCIENCES and partly from a survey conducted by the Con- 
necticut State Board of Fisheries and Game; the potentiometer used was 
purchased with a grant from the Sheffield Fund of Yale University. The 
generosity of these bodies in making our researches possible is here grate- 
fully acknowledged. Our best thanks are also due to Dr. O. A. Johannsen 
of Cornell University for examining our Chironomidae. 
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ON IDEALS OF DIFFERENTIAL POLYNOMIALS 
By J. F. Rit 
DEPARTMENT OF MATHEMATICS, COLUMBIA UNIVERSITY 
Communicated January 4, 1939 


Notwithstanding the excellent work of Raudenbush on perfect ideals, 
the general theory of ideals of differential polynomials is in a state which 
can hardly be called definitive. We hope to increase the plausibility of 
the existence of a complete theory in proving here a theorem of some 
generality on unrestricted ideals. Applied to the very special case of a 
single differential equation, in one unknown, of the first order, this theorem 
provides a substantial idealtheoretic background for the theory of the 
singular solutions of the equation. 


Ww 
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THEOREM: Let > be an ideal of differential polynomials in the unknowns 
Yi, ---» Yn» Let the manifold Mm of = be composed of manifolds, not neces- 
sarily irreducible, Mi, ..., Ms, none of which has a solution in common with 
any other. Then > can be, represented in one and in only one way as the 
intersection of s ideals 2, ..., 25, the manifold of >; being MN;. 

One sees that it is sufficient to treat the case of s = 2. As Mi and Me 
are mutually exclusive, there exists a relation 


A+B=1, (1) 


where A holds 9; and B holds M:. For some p, (AB)? isin =. We raise 
both sides of (1) to the 2pth power. In the first member of the resulting 
equation, let MM represent the sum of those terms in which the exponent 
of A exceeds p, and N the sum of the remaining terms. Then 4+ N= 1 
and MN is in 3. 

We shall prove that the derivative M' of M is in =. The derivative of 
MN, which equals M’(2N — 1), isin 2. Thus 2 contains M’(4N? — 1). 
Because N?M’ is in = contains M’. WN’ is also in 

Now let =; and 2 be, respectively, the ideals generated by 2 + M and by 
= + N. Then and have M, and Me, respectively, as manifolds. 
We shall prove that 2 1s the intersection of 21 and Yo. 

It suffices to show that if a form G is contained in 2, and in 2», G is con- 
tained in =. Such a G, because it is in 21, has, by what precedes, an ex- 
pression C+ DM, with AsGisin2:, DMisin Because M = 
1 — Nand Nisin LetD= E+ FNwithEin>. Because 
Y contains MN, > contains DM and hence G. 

We settle now the question of uniqueness. Let = be the intersection of 
ideals 2; and 33, whose manifolds are, respectively, Nt; and Mb. 

Let G be any form in 3}. Some power N‘ of N is contained in 3}. 
Then GN’ is in , hence in 3;._ Because N = 1 — M,Gisin 3;. Again, 
let H be any form in 3}. For some ¢, 5; contains M‘ and hence HM’. 
As M = 1 — Nand HN isin 3, =; contains H. We have proved that 3; 
and >; are identical. So also, then, are Sz and 33. 

What precedes shows that, for any value of s, each 2; is obtained by ad- 
joining a suitable form M; to = and taking the ideal generated by > + 
M,. The forms of 2; will all be of the type C + DM;with Cin >. One 
sees that, if 2 has a basis for which a single exponent can be used, the 
same will be true for 2;, and the least possible exponent for 2; will not 
exceed that for 2. 

The above work goes over, with no essential additions, to forms in 
several independent variables. 
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GENERALIZATIONS TO SPACE OF THE CAUCHY AND MORERA 
THEOREMS 


By MAXWELL READE AND E. F. BECKENBACH 
DEPARTMENT OF MATHEMATICS, RICE INSTITUTE 
Communicated December 15, 1938 
Introduction.—Consider the function 
w = f(s) = x(u, v) + ixe(u, v), 2 = u + iv, (1) 


defined and continuous in a simply connected domain D. A necessary and 
sufficient condition that f(z) be analytic in D is that the Cauchy-Riemann 
equations be satisfied there: 


Aw = 0, (2) 
where 

Os 


is a differential operator. From (2) we obtain 
2 
Ox)? = 0. (3) 
Jj = 


According to the Cauchy and Morera theorems, a necessary and suffi- 
cient condition that the continuous function f(z) in (1) be analytic in the 
simply connected domain D is that for each circle C lying in D, 


S fade = 0. (4) 
Now (4) may be considered to be an integral analogue of the differential 
condition (2); it implies 


cxj(u, v)dz]* = 0, (5) 


which is analogous to (3). 
The real functions 


xj = x;(u, 2), j = 1, 2,3, (6) 


defined and continuous in a simply connected domain D, will be said to 
define a surface S. A generalization of (3) to space is 


3 
0, (7) 
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which characterizes isothermic maps; that is, maps which are conformal 
except at points where E = G = 0. 
An analogous generalization of (5) to space is 


In this note we summarize some of the results which we have obtained 
in the study of equation (8). A more complete discussion will appear in a 
subsequent publication. 

We shall have frequent use for the following direct computation. Let 
the functions (6) have continuous partial derivatives of the mth order in 
a simply connected domain D, so that they admit finite Taylor expan- 
sions about each point (uo, v) of D, 


xj(u, v) = (cos 6 + sin | + o(r”),7 = 1,2,3, (9) 


oa Ou ov 


where the partial derivatives are evaluated at (uo, %), where o(r“) denotes 
a quantity (not always the same quantity) such that 


lim _ 


r—0 0, 


and where 
u— = 7 cos — % = rsin 0. 
Then, for the circle C with center (1%, vo) and radius 7 in D, the left-hand 


member of (8) assumes the form 


= 


fur bets a1 2” 


"CPC, ~ (2° ~ + * GO 


(2 +i 22-12) 


m k yk +2 


where 


and where the “C, are binomial coefficients, 


k! 
sl(k — s)! 
I. Characterization of Isothermic Maps.—THEOREM 1. [If the functions 
(6) have continuous partial derivatives of the first order in a simply connected 


domain D, then a necessary and sufficient condition that they map D iso- 
thermically on a surface S is that for each point (uo, vo) of D, 


= 
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3 
v)dz]? = o(r'), (11) 
where C is the circle with center at (uo, vo) and radius r in D. 

Proor: With m = 1 in (9), (10) becomes 


From (12), it follows that (7) is both necessary and sufficient for (11). 

II. Characterization of Those Isothermic Spherical Maps Which Do 
Not Map Circles on Circles——Lemma. If the functions (6) are not identi- 
cally constant, and if they map D isothermically on a surface S that lies on a 
sphere Q of finite radius, then a necessary and sufficient condition that they 
map circles on circles is that the first quadratic form of S have the expression 


c?(du? + dv?) 
[(w — uo)? + (v — + a*]?’ 


ds* = 


e>0,c> 0. (13) 


Necessity. Consider the stereographic projection of Q on the (x1, x2)- 
plane, and let the map of S be D’. The product of the transformation 
(6) and the stereographic projection maps D on D’, and carries circles in 
D into circlesin D’. It follows that the product of the two transformations 
is equivalent to a single linear transformation, whence a computation 
yields ds* in the form (13). 

SuFFIcieNcy. Let the stereographic projection of D on the sphere 
Q’ be S’, where Q’ is the sphere with center at (uo, 1, a — c/2a) and radius 
p = c/2a, and where the pole of projection is at (uM, v, a). The first 
quadratic form of S’ is found to be identical with that of S, so that S and 
S’ are applicable. Therefore S and S’ are either congruent or symmetric; 
in either case, since circles in D are mapped on circles on S’, it follows that 
circles in D are mapped on circles on S. 

We note that if one non-null circle in D is mapped on a circle on the 
spherical surface S by the isothermic functions (6), then all circles in D 
are mapped on circles on S. 

THEOREM 2. If the functions (6) have continuous partial derivatives of 
the third order in a simply connected domain D, then a necessary and sufficient 
condition that they map D isothermically on a surface S that les on a sphere 
of finite radius, such that circles are not mapped on circles, is that 


dsl? = of) (14) 


hold identically in D for a = 6, but not for a = 8. 


3 
j=1 
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NeEcsssity. With m = 3 in (9), (10) becomes 


3 3 


+ + 6%). (15) 
For our isothermic spherical map, we have 

E=G,F=0; e=g,f=0; e = k= const. #0, (16) 


where E, F, G and e, f, g are the coefficients of the first and second funda- 
mental quadratic forms for S. The formulas of Gauss, simplified by (16), 
yield 


3 
(Axj)(Adxj) = 0. (17) 
= 
It follows from (7), (17) and (15) that (14) holds for a = 6. 
Suppose that (14) holds identically in D for a = 8. Then, by (15), we 
have (7), (17) and 
3 3 
2 D0 + 3 (Adz)? = 0. 
j=1 j=l 
A computation, using (16), yields 
3 3 
j=l j=1 


whence 
3(AZE)? — 2E(\2E) = 0. (18) 
It follows that 
NE = (19) 


where = 0, so that 9(z) is analytic in D. From the imaginary part 
of the left-hand member of (18), we obtain 


3E,du 2Eywdu 3E,dv 2Eurdv 


whence 
NE = E*/2 (er) (20) 


It follows from (19), (20) and the Cauchy-Riemann equations, that 


= Apu + = + Qe, (21) 
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where do, @i, @2 are real constants. Substituting from (21) in (20), and 
integrating, we find 


4 = 
(u? + + ayu + aw + 


(22) 


Since on S the Gaussian curvature K is positive, it follows from 


1 
K=- OE A log E 
that (22) may be written in the form (13). Hence, by the above lemma 
circles in D are mapped on circles on S, contrary to hypothesis. 
Surricrency. Applying the differential operator \ to (7), we obtain 


3 
443) = 0, (23) 
Jj = 
and then applying \ to (23), and using (17), we obtain 
3 
= 0. (24) 
Jj = 


The four real linear homogeneous equations in Ax;, 7 = 1, 2, 3, implied 
by (23) and (24), have a solution other than Ax; = 0, 7 = 1, 2, 3, only 
if the rank ¢ of the matrix of coefficients of the system is less than three; 
but a computation shows that ¢ < 3 implies e = g, f = 0. Hence, (7) 
and (17) imply that D is mapped isothermically on either a minimal 
surface or a sphere. 

If D is mapped isothermically on a minimal surface, that is, if (7) and 
Ax; = 0,7 = 1, 2, 3, are satisfied, then, by (10), (8) is satisfied for each 
circle C in D; in particular, (14) is satisfied identically in D for « = 8, 
contrary to hypothesis. Or if D is mapped isothermically on a surface S 
that lies on a sphere Q, and if circles in D are mapped on circles on S, 
we consider any circle C in D and its image C’ on S, and project Q stereo- 
graphicaily on the plane II of C’, with D’ being the map of Don II; there 
is induced an isothermic mapping 


xj = yj(u,v), j = 1, 2, 3, 


of D on D’, that is, on a minimal surface; hence, 


ES cou, = 0; 


but on C, 


5 
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xj(u,v) = y,(u, 2), j = 1,2, 3, 
so that again (8) holds. 

Hence, under our hypotheses, D is mapped isothermically on a spherical 
surface of finite radius, and circles are not mapped on circles. 

III. Characterization of Those Isothermic Spherical Maps Which Map 
Circles on Circles and of Isothermic Maps on Minimal Surfaces —TueEo- 
REM 3. If the functions (6) have continuous partial derivatives of the third order 
in a simply connected domain D, then a necessary and sufficient condition 
that they map D isothermically either on a spherical surface, carrying circles 
into circles, or on a minimal surface, is that (8) hold for each circle C in D; 
further (14), with a = 8, is equivalent to (8). 

The proof is contained in the proof of Theorem 2. 


AN EQUIVALENCE THEOREM FOR SERIES’*OF ORTHOGONAL 
POLYNOMIALS 


By G. H. PEEBLES 


DEPARTMENT OF MATHEMATICS AND MECHANICS, INSTITUTE OF TECHNOLOGY, 
UNIVERSITY OF MINNESOTA 
Communicated January 9, 1939 
1. Introduction.—A set of polynomials [pn(x) = + ...] which 
satisfy the relations 
b 
_ j0,m n, 
p(X) Pn(x)Pm(x)dx = m =n, (1) 
b 
where p(x) is non-negative on (a, 0) and such that /p(x)dx > 0, may be 
a 


used to expand formally an ‘‘arbitrary’”’ function f(«): 


J b 


One may then consider the problem of showing for what interval, if any, 


n 
lim >> dybe(x) = f(x). 
k=0 

Concerning one of the several methods of attacking this problem, 
Shohat says,! ‘“‘The supreme goal, when dealing with [(2)], is to show that, 
regarding convergence, it behaves in a certain subinterval of (a, 0) like 
the ordinary Fourier series expansion of f(x) or of some function simply 
related to it. This constitutes what we call the equiconvergence theorem, 
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of greatest importance in our theory. Indeed, we well know the wide 
range of validity of the Fourier series expansion, which makes it such a 
powerful analytical instrument. The classes of [orthogonal polynomials] 
for which this goal is attainable are evidently the most interesting, and 
the method by which the goal has been attained is evidently the most 
powerful one.” 

This goal has been attained by A. Haar? and W. H. Young? for p(x) = 1 
(Legendre polynomials) and by G. Szegé‘ for a more general weight func- 
tion. The purpose of this paper is to present a compact and self-contained 
proof of equivalence, without the use of asymptotic formulas, which by 
reason of its simplicity seems to render the whole theory much more 
accessible, and at the same time, without reproducing all the results, of 
more elaborate investigations, leads incidentally to an extensive range of 
conclusions which are believed to be new. 

A polynomial expansion which is shown to be equivalent to a Legendre 
polynomial expansion is in turn equivalent to a Fourier series expansion, 
if the restrictions given by Haar and Young hold. A more direct con- 
nection, however, can be made with a Fourier series expansion by means 
of the so-called trigonometric polynomials. If the interval (a, 5) is taken 
for simplicity as (—1, 1), the trigonometric polynomials are those orthogo- 
nal with respect to the weight function (1 — x?)~'”. Relation (2) for 
these polynomials becomes by the substitution x = cos 6 the Fourier 
series expansion of the even function f(cos 6). 

The method used in proving the equivalence theorem will be seen to 
be applicable to systems of orthogonal functions other than polynomials.® 
It depends on a simple relation for the difference of two partial sums. 
Suppose the system of polynomials [p,(x; »:)] are orthogonal and normal- 
ized with respect to the weight function p;(x), the polynomials [pn(x; p2) | 
are orthogonal and normalized with respect to the weight function p2(x), 
where, here and hereafter, it is understood that p;(x) and pe(x) are non- 
negative, positive on a set of positive measure and of class L(a, b). Let 

n 


n 


Sn(X; pr) = Dipi(x; and sa(x; = aipi(x; p2) be the nth partial 
i=0 


sums of the formal expansions in terms of the two systems of a func- 
tion f(x) of classes Lp(a, b), L?p,(a, b) and Lp(a, 6). The difference 
Sn(X; p2) — Sn(x; pi) is a polynomial of degree n at most; so 


Su(x; p2) Su(x; p1) (Xx; p2); (3) 


k=0 


where rz, depends on n and is given by the formula 


b n 
= JS p2(x) [aipi(x; p2) — bips(x; pe)dx. 
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By virtue of (1), = Qe — b; pi(x; pi)Pe(x; p2)dx. Assume 
that a function p(x) of class Lp(a, °b) ota such that p(x)p;(x) = po(x) 
almost everywhere in (a, b). Then pro(x)pi(x; pilPe(x; = 
JS [p(x) Pe(x; pe) ]pi(x; pi)dx and, if Bir represents the last integral, one 
from Parseval’s formula = biBi = pi(x) p(x) pe)f(x)dx = ap. 


i=0 
Hence 


and lim fen = 0. 
Equation (3) with the relation just obtained enables one to place on 
p(x) additional conditions sufficient for lim [sn(x; p2) — sn(x; pi)] = 0 


wherever the polynomials [p,(x; p2)] are bounded. It follows that if one 
series converges to a value, diverges, is summable, the other must likewise 
converge to the same value, diverge, be summable; if one series is uni- 
formly convergent in an interval where the polynomials [P,(x; p2)] are 
uniformly bounded, the other is uniformly convergent in this interval. 
The powerful convergence tests for Fourier series expansions become, 
through the medium of the function f(cos @), instruments for determining 
the convergence of any polynomial expansion of f(x) which can be shown 
to be equivalent to the trigonometric polynomial expansion. The present 
treatment is believed to constitute, in particular, a valuable simplifica- 
tion in the theory of convergence and summability of Legendre series as 
well as of more general series of Jacobi polynomials. 

2. Notation.—Subsequent discussion will be simplified if certain 
abridgments are permitted. The symbol 2(x; \) will be understood to 
|e |x — ro|* x 
represent a function of the form —-.”. 

lb—a lb—a b-—a 


where always \i, ...,A4, 2 Oanda Sn <m<... <7, Sb, and the 
symbol tm(x; a polynomial (: ny (: ay (: of 
b-a b-—a b-a 
degreem = w+ ue+... + non-negative on (a,b). Ther’sare arbitrary 
points of the interval (a, b), but, in any one instance, if there is occasion to 
write several such functions, say, r(x; w(x; mm(x; the r’s are to be 
the same foreach function. It will be convenient, when setting down condi- 
tions preparatory to establishing some contemplated result, to regard (x; \) 


This has the advantage of per- 


as consisting of a single factor . ies 
—a 


J 
tm = 
t=n+1 
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mitting one to replace a set of relations such as > > 
hss eas by \> 2’ and, it will be seen, amounts to first proving the result 
for the simple function and then extending it to the more general function. 
Products such as r(x; \)a(x; d’) will be written in the form (x; \ + X’). 

Further simplification is effected by substituting the statement, “‘p2(x) 
is equivalent to p:(x),”’ for the statement, ‘‘the two expansions of f(x) in 
terms of the polynomials corresponding to the weight functions p2(x) and 
pi(x) are equivalent.” 

3. Preliminary Formulation.—A special form of relation (3) is required 
for the procedure given here. Let po(x) = (x; d)po(x) and without 
change in the restrictions on and definitions of p(x), pe(x), p(x) and f(x), 
but with slight modifications of the definitions of rx and Bjx, write relation 
(3) for pe(x) and p;(x): 


Sn(x; po) — Sa(x; pi) = TenPr(X; 


b 
Tin = di S d) p(x) po) pr)dx = 


t=n+1 
One sees from the definition of Bj that 
A)p(x) pe(x; p2) ~ Birpi(x; pi). (4) 


b 


The existence of / pi(x)[p(x)|°dx assures the existence of 
a 


b 
Hence 
b 
Bix pi(x) [w(x; d) p(x) Pe (x; pe) = 
b 
S a(x; d)p(x)p2(x) [De(x; p2)]*dx (5) 
and 


b 
> JS pr(x) [ar(x; A) p(x) P(x; p2) — on(x; k)]*dx, (6) 


t=n+1 


where on(x; k) = >> Bixpi(x; pi). 
i=0 


Let am(x; 4) be a polynomial of the kind described above, u < 2; and, 
for every value of k < n — m, let P,-~—~m(x) be a polynomial of degree 
n — k — m at most. Since the integral in equation (6) is increased, if 
on(x; k) is replaced by any other polynomial, one has, fork < n — m, 


n 
k=0 
n 


f 
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i=n 


6 
Bis s pi(x) [ar(x; A) p(x) pe (x; p2) — 
b 
S pr(x) — w)p(x) — (7) 
Suppose tm(x; and the P’s can be taken so that [am(x; u)]?pi(x) < 


w(x; A)po(x) and — w)p(x) — | S en—-k—m, a constant. 
Then 


i=n+1 
b 
S r(x; d)pe(x) [pe(x; = —e—m- 


By Cauchy’s inequality, 


f=n+1 j=n+1 j=n+1 


Ykn 


Hence, if M(x) = | pe(x; pe) | in (a, b) fork = 0, 1, ..., one has from (3) 


Su( Ba) — < aca) { > 


t=n+1 


+ Bis |}. (9) 


k=n—m+1Li=n+1 


There is some latitude in the choice of assumptions sufficient to cause 
the second member of the last inequality to tend to zero as m increases 
indefinitely. Since more satisfactory results are obtained, in general, if 
f(x) is lightly restricted, f(x) will be left, as heretofore, a function of classes 
Lp(a, b) and Lp(a, b) and the severest restrictions placed elsewhere. 

Suppose d/dx[r(x; — )p(x)] satisfies a Lipschitz condition of order 
aconstant,0 <a<1. Then, for m > 1, there exists a polynomial P,(x) 
of degree n at most such that | r(x; \ — u)p(x) — Pn(x) | < K/n'** in the 
interval (a, b), where K depends only on vy and the length of the interval.® 
Hence, the polynomials governing the magnitude of the e’s can be so chosen 


— K K K 
é) <a constant for all values of m. Equation (5) shows that Da B}, is 
t=n+1 


bounded for all values of m and k, if r(x; \)p(x) is bounded in (a, b). The 
number of such terms appearing in (9) is a fixed number m. The assump- 
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tion that f(x) is of class L*p,(a, 6) implies that lim = 0. There- 


t=n+1 


fore, if M(x») exists, lim [sn(xo; po) — Sn(x0; pi)] = 0. 


© 
As the boundedness of M(x) is important in this result, it is useful to 
know that the polynomials [p,(x; p2)] are bounded at a value of x which 
is not a zero of r(x; A)p(x) and for which the polynomials [f,(x; p:)] are 
bounded.’ The conditions given are sufficient to permit one to replace 
the sign of formal expansion in (4) by the sign of equality. Since By, = 


JS a(x; )po(x)pe(x; po) pi(x; pi)dx = 0, when i < k, w(x; d)p(x)pe(x; po) = 


© © 1/2 K 
2 Birpi(x; But | By | < {Esa} < when 


i—k—m-—1> 1. Hence, if N(x) = | pe(x; | fork = 0,1, ..., 


k+m+2 K K 
No(s) | pals; NG) Bal + + 


= 


ak ~ +.. : As Bix is uniformly bounded, the desired result 


follows immediately from this inequality. 

The conditions which are sufficient for these results are collected for 
future reference in a 

Lemma. [f f(x) is of classes Lpi(a, b) and L*p,(a, b) and tf there is a 
function p(x) and a polynomial r(x; such that in (a, b) = 
p(x)pi(x) almost everywhere, )p2(x) = [m(x; and d/dx[m(x; 
d — y)p(x)] satisfies a Lipschitz condition of order a, then the two expansions 
of £(x) in terms of the polynomials [p,(x; p2)| and [pa(x; p1)] are equivalent 
wherever the former system is bounded. The polynomials [pa(x; p2)] are 
bounded wherever the polynomials [p,(x; pi)] are bounded except for the zeros 
of d)p(x). 

4. Equivalence Theorem.—From the preceding Lemma it is possible 
to draw a succession of results which can be combined in a single theorem. 
The first result, however, comes more easily from the procedure leading 
tothe Lemma. Let r(x; \) = mm(x; u) = 1, and let p(x) be a polynomial of 
degree r, non-negative in (a,b). Thenfor m — k = ronecan take P,_,(x) = 
p(x) and so have from (7) ua Bi, = 0. The first inequality of (8) shows, 

7 1 


therefore, that rin = 0,k = 0,1,...,2 — 7, and, since }> B}, is bounded 


t=n+1 
for all values of and k, that rin, k =n —r + 1,...,, tends uniformly to 
zero as n increases indefinitely. One has from (3) p(x)pi(x) equivalent to 
pi(x) wherever the polynomials [pn(x; ppi)] are bounded. One also has 
the polynomials [p,(x; pp:)] bounded where the polynomials [p»(x; p1)] 
are bounded except for the roots of p(x). Hence, pi(x) and p(x)pi(x) are 
equivalent wherever the polynomials [p,(x; p:)] are bounded and p(x) ¥ 0. 


"oo 
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As a first application of the Lemma, let po(x) = pi(x) = (x; d’)pi(x), 
where p(x) has the properties hitherto attributed to p(x). In abridged 
notation, let the integer = \’ + land) = 2u —2d’. Thenz(x; d)po(x) = 
a(x; A’)pi(x) = + = w(x; Qu)pi(x) = w(x; Qu + 
= [rm(x; Since’ — = w— 1 and p(x) = 1, 
d/dx (r(x; % — u)p(x)] satisfies a Lipschitz condition of order a in (a, b). 
Hence, by the Lemma, p2(x) is equivalent to p:(x) wherever the polynomials 
[Pn(x; p2)] are bounded. But p(x) = [am(x; u)]%p:(x) is also equivalent to 
pi(x) wherever the polynomials [pn(x; p2)] are bounded. Therefore 
a(x; d’)pi(x) and pi(x) are equivalent wherever the polynomials [p_(x; p2)] 
are bounded, or, using part of the last result, wherever the polynomials 
[Pn(x; p1)] are bounded with the exception of the points 7, 72, ... 1,. 

To advance still further suppose po(x) = vm(x; d’)pi(x), v, a positive 
constant. In the Lemma let \ be an even integer greater than or equal to )’ 
and let = A. Then a(x; — u)p(x) = p(x) and pe(x) = (x; A)p2(x) = 
vm(x; + = ve(x; 2A)pi(x) = The Lemma calls 
for p2(x) = [rm (x; «)]%p:(x), whereas, here po(x) = v[am(x; u)]2:1(x). But 
this is readily seen to present no difficulty. Therefore, if p’(x) is assumed 
to satisfy a Lipschitz condition of order a, p2(x) and pi(x) are equivalent 
wherever the polynomials [pn»(x; p2)] are bounded. By the preceding 
result, the same statement can be made for po(x) and p(x). So pe(x) 
and pi(x) are equivalent where the polynomials [pn(x; :)] are bounded 
and r(x; ’)p(x) 0. 

The last case which includes most of the earlier results is contained in the 

TueoreM. If the polynomials [pa(x; p2)] and [pa(x; p:)] are orthogonal 
and normalized with respect to the weight functions p2(x) and p(x) over the 
interval (a, b), if £(x) is of classes Lp(a, b) and L*p;(a, b), and if a function 
p(x) and a set of points 1, to, ..., ¥,, exist such that in (a, b) p(x) = 


|x — 1; |™pi(x), po(x) = p(x)pi(x) almost everywhere, and p'(x) satis- 


fies a Lipschitz condition of order a, where v is a positive constant and the 
d’s are non-negative constants, then 


6 
lim p2) S pr(t)f(t)Pr(t; p2)dt — 


© 


r(x; pr) S pr(t)f(t)pe(t; pi)dt] = 0 


wherever the polynomials pi)] are bounded and p(x) x —1,|"0. 


This Theorem can be further generalized. If p(x) = (1 — x*)~¥?, 
p’(x) does not satisfy a Lipschitz condition. But the derivative of the 
product (1 — x?)”*o(x), after removing discontinuities, does. One has, 
then, (1 — x?)!/p.(x) equivalent to p:(x), (1 — equivalent to 
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p2(x) and, for a function of class Lpe(a, b), p2(x) equivalent to p:(x) wherever 
the polynomials [p,(x; p:)] are bounded except at —1 and 1. It is clear 
from this example that one may supplement the Theorem by the statement: 
The restriction on p'(x) may be replaced by two conditions, f(x) is of class 
Lp.(a, b) and a set of points rj, r3, ...., 1, exists such that in (a, b) d/dx 


[o(x) I |x — 1; |‘) satisfies a Lipschitz condition of order a, where the \’’s 
i=1 


are non-negative constants. The conclusion then becomes po(x) and p,(x) 
are equivalent wherever the polynomials [p,(x; pi)] are bounded and 


p(x) — |x #0. 
i=1 i=1 


1 Shohat, J., Bull. Am. Math. Soc., 41, 51 (1935). 

2 Haar, A., Math. Ann., 78, 121-136 (1918). 

3 Young, W. H., Comptes Rendus, 165, 696-699 (1917). 

4 Szegd, G., Math. Zeit., 12, 61-94 (1922). 

5 A general treatment, which sets forth necessary and sufficient conditions for equiva- 
lence of expansions in terms of orthogonal functions, is given by Walsh, J. L., and 
Wiener, N., Jour. Math. and Phys. of the Mass. Inst. Tech., 1, 103-122 (1922). 

6 For theorems on the order of approximation of continuous functions by polynomials, 
see Jackson, D., ‘“‘The Theory of Approximation,’”’ Am. Math. Soc. Colloquium, 11, 
Chap. I. 

7 The results on the boundedness of orthogonal polynomials to be obtained here are 
sufficient for the present needs. For a more complete, but less simple, treatment, see 
Bernstein, S., Jour. Math. pures et appliquées (9), 9, 127-177 (1930), and 10, 219-286 
(1931). 

8 See Jackson, D., Ann. Math., 34, 533 (1933), Theorem I. 

® Cf. Jackson, loc. cit.,8 p. 534, where further references are given. 
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I. Introduction.—We shall begin by giving some of the fundamental 
concepts in the geometry of curvature elements of the plane. By a 
curvature element, we mean simply a point of a curve together with the 
first two derivatives of the curve at the given point. A curvature element 
is given by the four numbers (x,.y, p, g), where (x, y) are the cartesian 
coérd‘nates of the point, p is the derivative of y with respect to x, that is, 
p = dy/dx, and q is the derivative of y with respect to x, that is, g = 


VoL. 25, 1939 MATHEMATICS: KASNER AND DE CICCO 105 


dp/dx = d*y/dx*. We observe that a curvature element may be repre- 
sented in the plane by a point together with a circular arc passing through 
that point. We shall call the totality of ©‘ curvature elements of the 
plane a plenum. ; 

A set of ©! curvature elements of a plane is called a series of curvature 
elements. A series may be given by the three equations y = y(x), p = 
p(x), q = q(x), where y, p, q are arbitrary functions of x only. If a series 
satisfies the two conditions dy—pdx = 0, dp-qdx = 0, the series is termed a 
union. A union obviously consists of a curve or a single fixed lineal 
element (x0, yo, fo). 

A collection of ©? curvature elements is said to be a field of curvature 
elements. A field may be given by the two equations p = p(x, v), gq = 
q(x, y), where p and gq are arbitrary functions of (x, y) only. If there 
exist ©! unions which fit the curvature elements of a given field (that is, 
all the curvature elements of the ©! unions coincide exactly with all the 
curvature elements of the given field), the field is said to be an integrable 
field. Of course not all fields are integrable fields. The necessary and 
sufficient condition that the field p = p (x, y), q = q (x, y), be an integrable 
field is that the functions p and q of (x, y) satisfy the single partial differential 
equation of first order q = Px + PpPpy- 

A set of 3 curvature elements is termed an opulence of curvature ele- 
ments. An opulence may be given by the single equation g = q(x, y, p), 
where q is an arbitrary function of (x, y, ») only. An opulence corre- 
sponds to a differential equation of the second order. We can always 
find ©? unions whose curvature elements coincide exactly with the curva- 
ture elements of a given opulence. Therefore every opulence is an in- 
tegrable opulence. 

A transformation between the curvature elements of two planes shall be 
termed a curvature element transformation. A curvature element trans- 
formation does not, in general, convert every union into a union. If a 
curvature element transformation converts every union into a union, then 
it is called a contact transformation. It is a well-known theorem of Back- 
lund that the group of contact transformations of curvature elements is iden- 
tical with the group of extended contact transformations of lineal elements. 

A curvature element transformation does not, in general, convert every 
integrable field of curvature elements into an integrable field of curvature 
elements. The problem of this paper is to find the group of curvature 
element transformations which convert every integrable field into an 
integrable field. Our fundamental result is that our group of transforma- 
tions is the group of contact transformations, that is, the group of extended 
contact transformations of lineal elements. Thus the solution of our problem 
gives us a new characteristic property of the group of extended contact 
transformations of lineal elements. 
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We note that we do mot demand that a single union of an integrable 
field shall be converted into a single union of the transformed integrable 
field. But it does result from our work that if a transformation carries 
every integrable field into an integrable field, then a single union of an 
integrable field is sent into a single union of the transformed integrable 
field. 

II. Integrable Fields of Curvature Elements into Integrable Fields of 
Curvature Elements —FUNDAMENTAL THEOREM. The necessary and suffi- 
cient conditions that the curvature element transformation 

X = X(x, 9, p,q), Y = 9), P = P(x, y, 9), Q = Q(x, 
converts every integrable field of curvature elements into an integrable field of 
curvature elements are that the functions X, Y, P, Q satisfy the equations 


x= X (x, p), Y= Y(x, p), 


tae, 


_ + pP, + 


+ ox, + aX,’ 


where X and Y are functions of (x, y, p) only. Thus our group of trans- 
formations is the group of contact transformations of curvature elements. In 
other words, it is the group of extended contact transformations of lineal elements. 
The sufficiency of our theorem is obvious. The remainder of the paper 
is concerned with the proof of the necessity of our theorem. 
Under our transformation, we must have 


whenever g = p, + Pp,. 


Q = Px + PPy, (1) 


The equation (1) may be written in the form 


Q X, + + Xy + Xppy + 
P, + Pyp. + Pa: Py + + Pay (2) 
Y,+ + Ya. Yy + Vopy + Yaa, 
4 Xe t+ + Xo Xy + + 
P, + Pop. + Pez Py + Pypy + Poy |’ 


where p, = q-pp,. The equation (2) must be an identity in py, g,, gy. 
Hence, upon setting the coefficients equal to zero, we obtain the five partial 
differential equations of first order 


= 


dy: 
ial 
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Q Y, Y, Ve Y, P, P, 
Xq Xx, P, P, | X, 
Y, Y, + P, P, : 
x, 
Y,+>Y, Y, Y, + pY, Y, P,+pP, P,|’ 
Q X, + xX, Py +P X, + xX, 
Y,+Y, Y, P,+qP, P, 


In the rest of the paper, we shall concern ourselves with the solution of 
these complicated equations. We shall show that the only possible solu- 
tion for the four functions X, Y, P, Q with non-vanishing jacobian is 
furnished by the equations (L), and therefore we have actually an ex- 
tended contact transformation. 

Now let a, b, c, d represent the quantities 


a = P, — QX,,b = Py — QXy,¢ = Py — QX,,d = P, — QX,. (4) 


Before proceeding, it must be observed that at least one of the four quan- 

tities a, b, c, d must be different from zero. For otherwise the jacobian of 

the transformation would be zero. This is contrary to our hypotheses. 
By means of (4), the five equations (3) may be put in the form 


a(Y, — PX,) — — PX,) = 0, 

a(Y, — PX,) — c(Y, — PX,) = 0, 

a(Y, — PX,) — d(Y, — PX,) = 0, (5) 
b[Y, + pY, — P(X, + pXy)] — (pc + d)(Vy — PX,) = 0, 

(qb + d)(Y, — PXy) — c[Y, + qY, — P(X, + = 0. 


We shall prove that Y, = PXq. For let Y,  PX,. Then we can solve 
the first three of equations (5) for }, c, d, obtaining 
(Y, — PX,)a (Y, — PX,)a wn (Y, — PX,)a 


Upon substituting the values of }, c, d, as given by (6), into the last two of 
equations (5), we find that they are identically satisfied. Hence under 
our assumption that Y, # PX,, we find that the equations (6) give the 
complete solution of the equations (5). 

From the above equations, we find that a ~ 0. For otherwise all four 
of the quantities a, b, c, d would be zero. Therefore by the assumption 


> 
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that Y, ¥ PX,, and by the fact that a ¥ 0, and by equations (4) and (6), 
we see that equations (6) are equivalent to the equations 
Py QX, — OX» Py-QXy Py — QX; (7) 
Y,— PX, Y,—PX,. Y,-PX, Y,-PxX, “ 


Thus from what was said above, the equations (7) are the complete solution 
of the equations (5) under the assumption that Y, + PX,. 

We observe that the jacobian J of our transformation can be written 
in the form 


x, x, x, x, 
P, Py Py P, OX, ~ OX, P, ~ OX, 
Q, Q, Q, 2, Q, QO, Q, Q, 
(8) 


By means of the equations (7), we find that the value of the jacobian J of 
the transformation, as given by (8), is zero. Thus we have proved our 
assertion that Y, = PX,. 

Since Y, = PX,, we find that the equations (5) become 


a(Y, — PX,) = 0,a(Y, — PX,) = 0,a(Y, — PX,) = 0, 
b[Y, + pY, — P(X, + pXy)] — (pbc + d)(Y, — PXy) = 0, (9) 
(gb + d)(Y, — PX,) — c[Y¥, + qV, — P(X, + qXp)] = 0. 


We shall show thata = Pg — QX, = 0. For if a # 0, then from what 
we have just proved and from the first three of the above equations, we 
obtain 

Y,= PX, Y, = PX,, ¥, = PXy, Y, = PX, (10) 


It is easily seen by these equations that the jacobian J of our transforma- 
tion vanishes. Since this is impossible, the validity of our assertion that 
a = 0 is established. 

Since we have proved that Y, = PX, and a = 0, the equations (5) 
reduce to the two equations 


DIY, + PY, — P(X, + PX,)] — +. PX) = 0 Gy 
(gb + d)(Y, PXy) AY, + — P(X, + »)] = 0. 


We now proceed to discuss the solution of (11). In order to simplify this 
discussion, we introduce the substitution 


h=d+ pe+ qb. (12) 


7 this 


12) 


VoL. 25, 1939 MATHEMATICS: KASNER AND DE CICCO 109 


Upon eliminating d from (11) and (12), we find that equations (11) are 
equivalent to the two equations 

WHY, + + oF, ~ FU, + My + qX»)] MY, — PX,) = 0, 

(13) 

clY, + pY, + q¥, — PX, + pX, + @X,)] — MY, — PX,) = 0. 

We shall prove that Y, + pYy + qYp — P(X, + pX, + qX,) = 0. For 
let us suppose that this is not the case. From the above equations, we 
find that \ ~ 0. For otherwise ) = c = 0, and hence from (12), 6 = ¢ = 
d = 0. Thus all four of the quantities a, 6, c, d are zero. This is im- 
possible. 

Now since Y, + pY, + qY, — P(X, + pXy + qX,) #0, 
we find upon dividing the two equations (13) and making use of (4) that 


Py Py — QXy 
¥,— Fx, 


By use of (4) and (12), we find from the first of equations (13) 


Ps t+ bPy + — Q(Xs + + (15) 
Y,— PX, Y,+p¥,+q¥, — P(X. + pXy + 


It follows from (14) and (15) that 


By these equations, the fact that Y, = PX, and the fact thata = P, — 
QX, = 0, we find that the jacobian (8) of the transformation is zero. This 


is impossible. Thus we have proved that 
Y, + pY, + = + pXy + 


Since Y, + pX, + qY, = P(X, + pX,y + qX,), the equations (13) then 
reduce to the equations 


MY, — PX,) = 0, — = 0. (17) 


We shall prove that} = P, + pPy + qP, — Q(X, + pXy + qX,) = 0. 
For if \ # 0, then from the fact that Y, = PX,, and from what we have 
just proved above and from the preceding equations (17), we obtain the 
equations (10). Thus the jacobian of our transformation vanishes so 
that it follows that \ = 0. 

We have thus far obtained the following solution for our system of 
equations (3) or (5) 


(14) 


Py = QXy + pPy + = + pXy + 
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The equations (18) obviously satisfy the equations (3) or (5) identically. 
It remains to prove that the only functions X, Y, P, Q with non-vanishing 
jacobian which satisfy (18) are given by the equations (L) of our theorem. 
We shall prove that X, = 0. Letus suppose that X, ~ 0. Then Y, +0 
and P, = 0. By eliminating P and Q from the equations (18), we find 


Y, + p¥y + = (Xz + + 4X5) 
(19) 


dx dy op) X, a dq \X,)° 
Upon taking the partial derivative with respect to q of the first of equa- 
tions (19), and also simplifying the second of equations (19), we obtain 
Ah PV + IV pq) V(X xq + 
PX + gX pq) + px, + dg 
° qd 
(20) 
+ P¥ + pq) — Vq(Xaq + PX yq + WX pq) = 


+ px, + qX X. 
q 
Upon subtracting these two equations, we find 


(21) 


Upon substituting (21) into (18), we find that under the assumption that 
X, # 0, the equations (18) are equivalent to 


Y, Y, 
P= a = ja Y, + pY, = P(X, + pX,), 
(22) 
P, 
Q = + + aPy = + PX, + 


where X, ¥ 0, Y, ¥ 0, P, ¥ 0. 


Under the existing conditions we see that X, # 0 and hence Y, ¥ 0. 
For if X, = 0, then Y, = 0 and hence X and Y are independent of p. 
From (22), we see that X, Y, P, and Q are independent of p. This makes 
the jacobian of our transformation zero. This is impossible so that we 
must have X, ¥ 0 and hence Y, + 0. 

From equations (22) and the fact that X, + 0 and Y, + 0, we find 


= 
qd 


at 


es 
ve 


VoL. 25, 1939 MATHEMATICS: KASNER AND DE CICCO 


Y, 
Y, + PY, = (X,+ pXy) 


(23) 


4 _ (Xz + bXy + 0 (V4 
Ps 1 3p) X X, dq \X, 
Upon taking the partial derivative with respect to of the first of equations 
(23), and also simplifying the second of equations (23), we obtain 


(XpV, — XyVp) + + PV yp) — + PX yy) = 


d /Y, 


(24) 


+ PYop) — Yo(Xap + PXyp) = + 0X) (32); 
Upon subtracting these two equations, we find 
(25) 


Substituting (25) into (22), we find that under the assumption that X, ¥ 0, 
the equations (18) are equivalent to 


(26) 


where X, ¥ 0, Pie ~ 0, Y, #0, Y, ¥0,P,# 0. From equations (26), it 
follows immediately that the jacobian of our transformation is zero. This 
proves that X, must be zero. 

Since X, = 0, it can be shown that the equations (18) are equivalent to 
the equations (L) of our theorem. Thus our curvature element trans- 
formation must be an extended contact transformation of lineal elements. 
This completes the proof of our theorem. 

An analogous problem in space dealing with the transformation of lineal 
elements of curves and surfaces is discussed in another paper.! 

1 See Kasner, ‘Lineal Element Transformations Which Convert Normal Families of 


Curves into Normal Families,”” Duke Mathematical Journal, March (1939). The result 
gives a new characterization of space contact transformations. 
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